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ABSTRACT 


The  question  of  the  possible  existence  of  the  quantum  bound  states  (localized  states)  due  to 
special  geometries  has  been  a  long  standing  problem  in  quantum  theory.  Knowledge  of  quantum 
bound  states  in  crossed  nanowire  system  is  very  important  in  understanding  the  properties  of 
spatially  confined  nanostructures  which  are  promising  candidates  for  device  applications  such  as 
transistors,  amplifiers,  switches,  biosensors,  photo-detectors,  solar  cells,  lasers  and  light-emitting 
diodes.  This  study  focuses  mainly  on  investigating  the  angular  dependence  of  the  lowest  bound 
energy  state  for  the  model  system  of  an  electron  trapped  at  the  intersection  of  two  identical 
narrow  channels  (nanowires)  crossed  at  an  arbitrary  angle.  When  the  channels  are  perpendicular, 
such  a  classically  unbound  system  is  known  to  possess  a  quantum  bound  state.  We  used  the 
variational  principle  to  obtain  the  general  criterion  to  study  the  role  of  tilted  geometry  for  the 
existence  of  bound  state  of  an  electron  in  such  a  quantum  system.  Using  suitable  trial  wave 
functions,  calculations  were  carried  out  to  estimate  the  upper  bound  for  the  ground  state  energy 
of  an  electron  located  inside  the  crossed  nanowire  system.  The  results  of  our  calculations  show 
that  the  bound  state  energy  varies  as  the  squared  sine  of  the  intersection  angle  of  the  crossed 
nanowires.  We  have  found  that  the  system  is  strongly  bound  when  the  intersection  angle  is  90° 
and  the  system  is  unbound  when  the  intersection  angle  is  0°.  These  particular  features  of  the 
crossed  nanowire  system  can  be  exploited  to  design  single  electron  ultra-sensitive  switching 
devices  using  newly  developed  lithographic  and  etching  techniques  in  nanotechnology. 
Furthermore,  our  model  may  be  useful  to  interpret  electron  transport  peculiarities  in  realistic 
systems  such  as  semiconductor  nanowire  films,  nanorods  and  carbon  nanotube  bundles. 
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I.  INTRODUCTION 

A.  Literature  Review:  History  and  Background 

The  current  progress  in  lithographic  and  etching  techniques  in  semiconductor  technology 
make  it  much  easier  to  design  low  dimensional  quantum  systems  such  as  quantum  wells 
(two  dimensions),  quantum  wires  (one  dimension)  and  quantum  dots  (zero  dimension).  The 
smallness  yet  finiteness  of  these  structures  invoke  quantum  mechanical  effects  into  play.  The 
motion  of  an  electron  in  a  narrow  crossed  channels  has  been  studied  previously  by  Schult, 
Ravenhall  and  Wyld  [1].  This  paper  report  the  detailed  calculations  designed  to  determine 
whether  quantum  effects  can  cause  the  trapping  of  single  electron  at  the  intersection  of 
wires.  The  potential  shape  have  been  ignored  but  considered  the  channels  such  that  the 
potential  is  zero  within  the  cross  region  and  infinite  outside.  Schult  et  al.  solved  the  free 
particle  Schrodinger  equation  using  two  methods  for  the  square  geometry  [Fig.  1] . 


FIG.  1.  The  geometry  of  2D  crossed  wires  of  same  width  (w)  intersecting  at  right  angle. 

They  have  computed  the  energy  and  the  wave  function  for  an  electron  caught  at  the 
intersection  of  the  two  narrow  channels  intersecting  at  right  angle.  Their  first  model  was  a 
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FIG.  2.  Properties  of  a  crossed-wire  system,  as  a  function  of  kw,  with  w  being  the  width  of  the 
wire  and  the  energy  given  by  E  =  The  bound  states  are  shown  as  vertical  lines.  The  first 

vertical  line  show  the  bound  state  energies  E\  and  the  second  one  represent  excited  states.  The 
dotted  vertical  lines  are  for  the  propagation  threshold  energies  for  the  ground  and  the  excited 
states  respectively.  The  functions  are  R,  the  reflection  probability  in  the  incident  wire;  T,  the 
straight-through  transmission  probability;  S,  the  probability  of  sideways  transmission  into  the  side 
wires;  and  S' ,  the  probability  of  sideways  scattering  into  the  second  transverse  band,  n  =  2.  Figure 
is  courtesy  of  Schult  et  al.  [1] . 


mesh  point  method  in  which  they  replaced  the  2D  Schrodinger  equation 

Hty(x,y)  =  E^(x,y) 


(1) 


by  a  difference  equation  for  the  wave  function  evaluated  on  a  rectangular  mesh  of  points 
in  the  x  —  y  plane  with  discrete  evolution  in  pseudo-time  variable,  t.  They  treated  the 
Schrodinger  equation  as  the  discretization  of  the  differential  equation 


It 


=  V2T. 


(2) 


Taking  the  initial  wave  function  to  be  symmetric  about  the  center  lines  in  the  x  and  y 
directions  they  computed  the  energy  value  to  be  Ex  =  0.66 Et.  Et  is  the  threshold  energy 
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FIG.  3.  The  bound  state  wave  functions  for  the  ground  state.  The  wave  function  has  steep  slope 
near  the  corner,  concave  upward  along  the  channel  and  convex  upwards  across  the  channel.  Figure 
is  courtesy  of  Schult  et  al.  [1] . 

for  free  propagation  in  the  open  channel,  which  is  given  by 

=  ttHi2 

*  2mw 2  (  ’ 

where  w  is  the  width  of  the  channel  and  m  is  the  effective  mass  of  an  electron.  Also  for 
the  initial  wave  function  which  is  odd  about  the  x  and  y  center  lines  and  even  about  the 
x  =  y  line,  the  next  bound  energy  level  is  at  E-2  =  3.72 Et,  the  propagation  threshold  being  at 
kw  =  tv  [Fig.2],  This  mesh  point  method  gave  only  the  existence  of  two  bound  states.  Schult 
et  al.  used  eigenfunctions  expansion  to  solve  the  differential  equation  to  get  higher  order 
bound  states.  The  bound  states  in  case  of  bent  wires  in  the  form  of  ”L”  and  ”T”  shaped 
geometry  have  been  observed.  They  found  that  the  wave  functions  near  the  corners  of  the 
intersection  varies  as  p2/3,  where  p  is  the  distance  from  the  corner  [Fig. 3].  The  presence 
of  corners  formed  by  the  intersection  of  two  nanowires  is  responsible  for  trapping  of  the 
electron.  The  electron  is  localized  in  the  square  region  and  cannot  propagate  along  the 
channel.  Furthermore,  they  estimated  the  Coulomb  energy  for  the  crossed  wires.  Using  the 
unperturbed  wave  function,  they  showed  that  two  electron  cannot  be  bound  in  the  cross 
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wires  because  Coulomb  repulsion  is  much  greater  than  (about  10  times)  the  single-particle 
binding  energy.  One  of  the  electron  can  be  bound  while  the  second  one,  even  with  opposite 
spin  is  pushed  into  the  subband  continuum. 


FIG.  4.  The  geometry  of  the  crossbar  structure  connecting  left  and  right  reservoirs  with  2D 
electrons  with  equal  Fermi  energy.  The  crossbar  appear  to  an  electron  at  the  bound  state  energies 
nearly  infinite. 

Berggren  and  Ji  tried  to  explore  the  possibility  of  observing,  at  least  in  principle,  the 
bound  states  in  the  two  intersecting  crossed  wires  by  resonant  tunneling  [2],  The  bound 
states  revealed  themselves  in  the  electric  transport  as  resonances  or  antiresonances  and  the 
tunneling  in  the  two  intersecting  wires  occur  through  discrete  states.  The  geometry  that 
they  considered  consisted  of  a  channel  connected  to  two  reservoirs  serving  as  a  source  and 
drain.  The  vertical  channel  was  replaced  by  a  bar  of  finite  length  [Fig.4],  They  suggested 
that  measurements  of  the  tunneling  current  could  be  one  of  the  ways  of  observing  bound 
states  in  the  cross  wires  .  For  such  a  crossed  region,  the  conductance  have  been  found  to  be 
quantized  as 

G  =  N(2e2/h)  (4) 

where  IV  =  1  or  2.  The  variation  of  the  conductance  as  a  function  of  Fermi  energy  show 
sharp  resonances  at  E'  and  E"  which  corresponds  to  the  localized  states  [Fig. 5].  The  two 
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FIG.  5.  The  conductance  as  a  function  of  the  Fermi  energy.  The  length  of  the  channel  connecting 
2D  electrons  and  its  width  are  kept  constant.  The  sharp  peaks  E'  and  E"  show  bound  state 
energies.  Figure  is  courtesy  of  Berggren  et  al.  [2]. 

peaks  corresponds  to  the  single  electron  resonant  tunneling  via  bound  state  energies  E\  and 
E2  as  estimated  by  Schult  et  al.  [1]. 

Soon  after  the  results  by  Schult  et  al.  for  the  one  electron  bound  states  in  an  intersecting 
crossed  wires,  Berggren  and  Wang  studied  the  possibility  of  two  electrons  being  trapped 
in  a  finite  cross  wires  intersection  [3].  The  structure  that  they  considered  consisted  of 
an  infinite  horizontal  channel  intersected  by  a  finite  perpendicular  channel  [Fig.6].  Their 
method  involved  was  solving  2D  Schrodinger  equation  self-consistently  for  each  electron 
taking  into  account  of  the  Coulomb  potential  of  the  other  assuming  helium  atom  model. 
They  found  that  the  first  electron  in  the  cross  region  experiences  an  effectively  attractive 
Coulomb  potential  while  the  second  is  expelled  from  the  central  region.  This  means  that 
no  self  consistent  solution  exist  for  the  second  electron  providing  the  evidence  of  Coulomb 
blockade  responsible  for  the  scattering  state  of  two  electrons  [Fig.  7].  Their  work  supplement 
the  earlier  result  obtained  by  Schult  et  al.  [1], 

Looking  back  to  the  previous  work  on  the  transport  properties  of  crossed  wires,  Ji  and 
Berggren  observed  that  the  quantum  bound  states  in  the  2D  semiconductor  structure  in 


FIG.  6.  (a)  Schematic  view  of  two  intersecting  quantum  wires  with  hard  walls.  The  width  of 
the  wires  is  w.  (b)  View  of  the  modified  quantum  cross  consisting  of  two  horizontal  leads  and  a 
perpendicular  bar  with  hard  walls. 


FIG.  7.  The  Coulomb  potential  seen  by  the  two  electrons:  (a)  The  Coulomb  potential  for  the  first 
electron  is  taken  as  attractive  potential.  In  (b),  the  Coulomb  potential  for  the  second  electron  acts 
as  a  potential  barrier. 

the  form  of  two  crossed  channels  intersecting  at  right  angle  give  rise  to  resonant  tunneling 
[4],  The  channels  that  they  considered  were  finite  and  one  of  them  connect  the  reservoirs 
of  2D  electron  gas  on  the  left  and  right  [Fig. 8].  The  potential  difference  was  applied  across 
the  two  reservoirs  and  conductance  values  were  measured  in  a  structure  with  the  number  of 
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FIG.  8.  The  geometry  of  the  number  of  crossbar  structures  connecting  left  and  right  reservoirs 
with  2D  electrons  with  equal  Fermi  energy.  When  the  potential  difference  is  applied,  electrons 
flow  from  left  to  right.  The  dashed  lines  represent  the  boundaries  where  the  wave  functions  are 
matched. 

intersections.  For  the  single  intersection,  the  conductance  showed  a  sharp  peak  below  the 
lowest  subband  threshold  Et.  This  peak  corresponds  to  the  tunneling  through  the  localized 
bound  states  at  the  center  of  the  cross  in  case  of  two  infinite  intersecting  channels  [1] .  They 
extended  the  calculation  to  a  double  cross  so  that  the  tunneling  below  the  threshold  can  be 
assumed  to  be  that  through  two  quantum  dots,  one  held  at  each  intersection.  These  quantum 
dots  are  coupled  by  a  weak  potential  and  the  resonance  split  into  two  peaks  given  by  Ei  and 
E2  [Fig. 9].  In  the  case  of  number  of  crosses,  the  resonant  tunneling  can  be  approximated 
by  Kronig-Penny  model.  There  are  as  many  resonance  peaks  below  the  lowest  subband 
threshold  as  there  are  crosses  in  the  structure  connecting  two  2 D  electron  reservoirs.  The 
calculated  value  of  the  conductance  reach  the  maximum  value  when  the  quantum  dots  are 
well  localized  to  an  intersection. 

So  far,  the  geometries  considered  were  2 D  crossed  wires  intersection  enclosing  square 
region  with  both  arms  either  infinite  [1],  semi-infinite  [3,  4],  or  finite  in  length  [2],  We 
have  seen  the  role  of  bound  states  in  the  electronic  transport  properties  of  two  dimensional 
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FIG.  9.  The  calculated  conductance  at  zero  temperature  as  a  function  of  Fermi  energy  E.  E\,  E% 
and  E3  denote  the  energies  at  which  resonant  tunneling  via  states  resembling  quantum  dots  takes 
place.  Figure  is  courtesy  of  Ji  et  al.  [4], 


crossed  structures  and  some  of  the  methods  of  observing  the  localized  states.  No  external 
held,  such  as  electric  or  magnetic,  was  applied  to  the  crossed  wires  system. 

Ravenhall,  Wyld  and  Schult  took  another  approach  to  probe  further  into  the  same  prob¬ 
lem  but  in  the  presence  of  external  magnetic  held.  Hence,  they  studied  rigorously  the  effect 
of  magnetic  held  on  the  bound  states  at  the  crossed  junction  of  nanowires  under  ballistic 
approximation  in  which  they  neglected  the  electron-electron  interaction  [5].  The  crossed 
wires  intersect  at  right  angle  and  the  magnetic  held  was  applied  perpendicular  to  the  plane 
of  the  intersecting  wires  [Fig. 10].  The  method  they  used  was  an  extension  for  the  zero  held 
case  to  the  case  of  a  magnetic  held  of  arbitrary  strength.  In  their  model,  the  voltage  that 
appear  across  top  and  bottom  of  the  wire  along  y-direction  gave  the  Hall  voltage  (Vyy).They 
calculated  the  current  (Ixx)  along  x-direction  which  depends  strongly  on  the  bound  states 
at  the  junction.  The  Hall  resistance  was  obtained  using 
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FIG.  10.  The  four-terminal  junction  used  to  study  the  effect  of  magnetic  fields  on  the  bound  state 
at  the  intersect. 

and  compared  with 

Rh  =  —2 
n0e2 

where  no  is  an  integer.  The  ballistic  approximation  predicts,  even  for  zero  field,  bound 
states  at  the  junction  and  also  the  variation  of  the  scattering  probabilities  with  energy  at 
the  threshold  for  the  opening  of  new  channels.  They  observed  that  with  non-zero  magnetic 
field  new  structures  are  associated  with  virtual  and  resonant  states  located  at  the  junction 
that  becomes  2 D  Landau  levels  for  large  value  of  the  magnetic  field  [Fig.  11], 

Immediately  after  the  results  for  2 D  system  of  crossed  wires  were  published,  the  quantum 
transmission  properties  of  an  electron  in  a  crossed  junction  of  two  out-of-plane  wires  have 
been  numerically  investigated  by  Takagaki  and  Ploog  [6].  They  considered  3 D  junction  with 
two  identical  rectangular  cross  wires  connected  at  right  angles  on  top  of  each  other  [Fig. 12], 
They  assumed  zero  potential  within  the  channels  and  infinite  outside,  so  that  for  a  given 
mode  (i,j),  the  threshold  for  the  propagation  is  given  by  (i2  +  j2)Et.  In  such  a  structure 
of  crossed  junctions,  an  electron  can  be  transmitted  in  the  direction  perpendicular  to  its 
incidence.  Their  calculations  show  that  when  these  wires  are  not  coupled,  the  electron  gets 
simply  transmitted  in  the  forward  direction  without  any  reflection.  However,  if  the  coupling 
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kw 

FIG.  11.  The  dimensionless  Hall  resistance  as  a  function  of  kw  for  dimensionless  magnetic  field 
B  =  6.0.  Figure  is  courtesy  of  Ravenhall  et  al.  [5]. 

is  non- zero,  the  electron  gets  reflected  with  probability  (Tr)  or  transmitted  in  the  forward 
direction  with  probability  (Tp)  and  into  each  end  of  the  upper  wire  with  probability  (Tg) 
such  that 

Tf  +  2  Ts  +  Tr  =  N ,  (7) 

where  N  «  nEp/4:Et  is  the  number  of  propagating  modes  in  the  channels  for  an  electron 
moving  with  Fermi  energy  EF.  The  transmission  exhibits  resonance  due  to  quasibound 
states  in  the  intersection  region  even  if  the  device  consist  of  straight  wires.  In  the  2D  in¬ 
plane  crossed  wires  junction  the  incident  mode  does  not  couple  with  the  bound  states  due 
to  parity  miss-match  hence,  no  resonances  occur  in  the  transmission.  The  appearance  of 
transmission  resonances  due  to  quasibound  states  gave  the  strong  evidence  of  the  quantum- 
mechanical  behavior  of  the  scattering  phenomena  (see  Fig. 13).  Furthermore,  they  calculated 
the  bend  resistance  Rr,  which  is  always  negative  for  2D  crossed  junction,  can  have  positive 
value  around  the  resonances  in  the  3 D  junction  (see  Fig.  14). 

An  idealized  geometry  of  3 D  junction  with  a  single  crossing  studied  by  Takagaki  and 
Ploog  have  limitations  with  regard  to  the  realistic  system  associated  with  number  of  cross¬ 
ings.  Recent  experiments  on  carbon  nanotube  crossed  with  a  semiconducting  one  have 
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FIG.  12.  Schematic  view  of  the  out-of-plane  crossed  wires  junction.  Two  rectangular  waveguides 
are  connected  at  right  angles  on  top  of  each  other.  The  width  of  the  waveguides  are  assumed  to 
be  equivalent.  An  electron  is  injected  from  the  lower  waveguides. 

shown  the  existence  of  bound  states  at  the  crossing  which  are  not  due  to  disorder.  With 
an  objective  at  clarifying  this  problem,  Makogon,  Jeu  and  Smith  made  a  detailed  study  of 
tunneling  effects  in  crossed  ID  systems  in  the  presence  of  potential  barriers  for  the  massive 
quasiparticle  excitations  [7].  Makagon  et  al.  considered  a  system  composed  of  two  layers  of 
crossed  quantum  wires  with  interlayer  coupling  separated  by  a  distance  (see  Fig. 15). 

Apart  from  ballistic  approximation,  they  took  into  account  of  the  electron-electron  inter¬ 
actions  in  their  calculations  to  study  the  tunneling  properties  in  different  cases  such  as  -  free 
electron  case,  single  crossing  of  two  crossed  wires  and  more  general  case  that  include  inho¬ 
mogeneous  potential.  In  the  case  of  single  crossing  of  two  crossed  wires,  they  observed  that 
in  addition  to  the  scattering  states,  bound  state  solutions  were  possible.  The  appearance 
of  bound  states  was  due  to  the  presence  of  tunneling.  However,  they  obtained  the  bound 
states  energy  much  smaller  than  computed  numerically  by  Schult  et  al. 

As  reported  earlier  in  a  classic  paper  by  Schult  et  ah,  corners  in  the  crossed  wire  system 
acted  as  traps  which  were  responsible  for  the  existence  of  localized  states  [1].  This  feature 
was  further  studied  by  Avishai,  Bessis,  Giraud  and  Mantica  in  the  case  of  open  geometries 
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FIG.  13.  Transmission  probabilities  as  a  function  of  energy.  The  thin  lines  represent  the  classical 
transmission  probabilities.  Figure  is  courtesy  of  Takagaki  and  Ploog  [6]. 


FIG.  14.  Bend  resistance  as  a  function  of  energy.  The  classical  of  Rb  with  the  zero-point  motion 
correction  is  represented  by  dotted  line.  The  arrow  indicate  the  propagation  threshold  for  modes 
The  inset  illustrate  the  current  and  voltage  probe  configuration.  Figure  is  courtesy  of 
Takagaki  and  Ploog  [6]. 
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X,  x2  x3 


FIG.  15.  Two  dimensional  array  of  crossed  wires  separated  by  a  distance  d. 


FIG.  16.  Geometry  used  for  the  existence  of  a  bound  state. 

[8].  In  this  paper  Avishai  et  al.  reported  the  problem  of  quantum  bound  states  in  open 
system  with  geometrical  constraints  on  free  propagation  without  any  explicit  potentials. 
In  such  case,  the  Schrodinger  equation  was  reduced  to  a  pure  Helmholtz  equation  with 
just  the  kinetic  energy  operator.  The  constraints  imposed  by  the  boundary  conditions 
on  the  wave  functions  were  strong  enough  to  generate  bound  states  in  a  completely  open 
geometries  which,  classically,  do  not  show  any  trapped  states.  They  considered  two  2 D 
semi-infinite  wires  joined  by  a  sharp  angle  known  as  ”  broken  strips”  (see  Fig.  16).  They 


16 


FIG.  17.  The  linearity  of  (7 r2  —  E)1!2  as  a  function  of  tan2  e.  Figure  is  courtesy  of  Avishai  et  al. 
[8], 

used  the  variational  principle  to  calculate  the  energy  of  an  electron  in  the  bent  strip.  Their 
result  proved  the  existence  of  at  least  one  bound  state  for  an  arbitrarily  small  angle  and  the 
binding  energy  was  found  to  be  proportional  to  the  fourth  power  of  the  sine  of  the  bent  angle. 
They  extended  their  geometry  for  the  case  of  a  crossed  wires  at  right  angle  and  calculated 
the  energy  (see  Fig.  17).  They  estimated  numerically,  the  upper  bound  to  the  ground  state 
energy  as  E  =  0.947 r2  as  compared  to  the  energy  E  =  0.927T2  obtained  by  Schult  et  al.  and 
0.937T2  by  Exner  and  Seba  [9]. 

In  order  to  apply  the  phenomenon  of  quantum  bound  states,  Takagaki  and  Ploog  consid¬ 
ered  the  more  realistic  situation  through  atomistic  approach  to  the  crossed  junction  problem 
[10].  They  investigated  the  quantum  transport  properties  in  cross  junctions  delineated  in 
an  anisotropic  two  dimensional  electron  gas.  They  modeled  their  anisotropic  system  using 
square  tight-binding  lattice  with  nonidentical  nearest-neighbor  hopping  amplitudes  in  the  di¬ 
rection  that  are  orthogonal  to  each  other.  The  model  they  used  consisted  of  a  symmetrically 
shaped  crossed  junction  defined  by  a  square  lattice  inclined  with  respect  to  the  anisotropy 
direction  (see  figure  Fig. 18).  They  assumed  that  the  probabilities  of  an  electron  turning 
into  the  left-hand-side  and  right-hand-side  leads  are  different.  Hence,  they  calculated  the 
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FIG.  18.  Geometry  for  tight-binding  model  for  anisotropic  cross  junctions  when  the  number  of 
the  transverse  lattice  sites  N  =  2.  The  solid  thick  lines  indicate  boundary  of  the  cross  junction. 
Takagaki  and  Ploog  [10]. 


FIG.  19.  Energy  dependence  of  transmission  probabilities  in  a  symmetrically  shaped  cross  junctions 
in  the  absence  of  magnetic  held.  Figure  is  courtesy  of  Takagaki  and  Ploog  [10]. 
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FIG.  20.  (a)  The  original  Y  junction  of  Q1D  channels  considered  here,  (b)the  reformulated  Y 
junction  of  tight-binding  channels,  and  (c)  the  effective  Y  junction  of  ID  channels  affiliated  with 
connection  scheme. 

transmission  probabilities  quantum-mechanically  using  tight-binding  model  and  found  that 
transmission  into  the  right-hand-side  lead  was  larger  than  that  into  the  left-hand-side  lead 
for  lower-lying  modes.  Takagaki  et  al.  observed  that  anisotropy  played  the  role  of  breaking 
the  symmetry  of  the  electronic  state  and  the  transport  symmetry,  giving  rise  to  a  transmis¬ 
sion  resonance  through  a  quasibound  state  in  the  cross  junction  and  a  finite  Hall  resistance 
despite  the  absence  of  magnetic  held  respectively  (see  Fig. 19).  The  magnitude  of  Hall  re¬ 
sistance  was  found  to  depend  weakly  on  the  number  of  occupied  subbands  in  the  leads. 
The  most  essential  features  of  Hall  resistance  they  observed  was  switching  of  the  polarity 
of  the  Hall  resistance  associated  with  the  subband  thresholds  that  occurred  in  the  quantum 
mechanical  regime. 

The  three-leg  junction  of  one  dimensional  channels  (Y  junction)  have  been  studied  by 
Voo  et  al.  (see  Fig.20).  In  such  geometry,  they  proposed  a  scheme  to  connect  the  wave 
function  on  different  one-dimensional  (ID)  branches  of  the  junction  [11],  Their  model  took 
into  account  of  the  difference  in  the  widths  of  the  quasi-one-dimensional  (Q1D)  channels 
in  different  systems.  They  tested  their  scheme  by  comparing  results  from  a  doubly  con- 
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FIG.  21.  The  transmission  probability  T  is  plotted  versus  the  dimensionless  longitudinal  wave 
number  in  the  case  of  quasi-one-dimensional  channels  and  in  case  one  dimensional  channels  for 
identical  channel  widths.  Figure  is  courtesy  of  Voo  et  al.  [11]. 

nected  one-dimensional  system  and  a  related  quasi-one-dimensional  system.  The  scheme 
they  suggested  may  be  useful  in  constructing  one-dimensional  effective  models  out  of  quasi- 
one-dimensional  channels.  Voo  et  al.  used  tight-binding  model  to  reformulate  the  two- 
dimensional  Schrodinger  equation  and  imposed  the  condition  of  continuity  of  wave  function 
at  the  N-leg  junction  of  one  dimensional  channels.  They  calculated  the  transmission  prob¬ 
abilities  for  the  Q1D  channels  and  ID  channels  as  a  function  of  the  wave  number  for  the 
case  identical  legs  of  the  Y  junction  (see  Fig. 21).  The  two  results  obtained  from  their  model 
agreed  with  each  other,  although  the  strength  of  the  attractive  nature  of  the  Y  junction  was 
underestimated. 

In  summary,  the  question  of  the  possible  existence  of  quantum  bound  (trapped,  local¬ 
ized)  states  by  special  geometry  has  been  a  long  standing  problem  in  quantum  theory.  In 
particular,  different  geometries  of  the  crossed  nanowire  intersections  explored  earlier  both 
in  two  and  three  dimensions,  have  shown  the  existence  of  bound  states.  The  ways  of  ob¬ 
serving  these  states  by  means  of  measurement  of  conductance  were  proposed  theoretically 
using  two  dimensional  electron  gas  models.  The  effects  of  the  external  magnetic  field  on  the 
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cross-wire  system  were  also  studied  which  helps  to  develop  scheme  of  observing  the  localized 
states  through  the  measurement  of  the  quantum  Hall  resistance.  The  presence  of  the  bound 
states  at  the  crossed  nanowire  intersections  greatly  affect  the  electron  tunneling  properties. 
Also,  quantum  bound  states  are  responsible  for  observed  anomaly  in  the  transmission  and 
reflection  properties  in  crossed  nanowires  systems  such  as  carbon  nanotube  bundles  and 
films.  Many  literature  sources,  till  now,  looked  into  the  effect  of  localized  states  on  the  elec¬ 
tron  tunneling  behavior  for  nanowires  crossed  at  right  angle,  but  failed  to  take  into  account 
the  dependence  of  the  localized  state  binding  energy  on  the  ( arbitrary )  intersection  angle. 
Hence,  the  objective  of  this  thesis  work  is  to  analyze  the  angular  dependence  of  the  bound 
state  energy  for  the  model  system  of  an  electron  trapped  at  the  intersection  of  two  identical 
narrow  channels  (quantum  wires)  crossed  at  an  arbitrary  angle  [12].  When  the  channels 
are  perpendicular,  such  a  classically  unbound  system  is  known  to  possess  a  quantum  bound 
state.  In  what  follows,  we  set  up  the  model  for  the  problem  and  investigate  the  angular 
dependence  of  the  electron  bound  state  energy  in  the  oblique  cross-wire  system. 
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II.  STATEMENT  OF  THE  PROBLEM 

A.  Objective 

The  main  objective  of  the  present  work  is  to  analyze  the  angular  dependence  of  the  bound 
state  energy  for  the  model  system  of  an  electron  trapped  at  the  intersection  of  two  identical 
narrow  channels  (quantum  wires)  crossed  at  an  arbitrary  angle  (see  Fig.22).  The  potential 
is  zero  inside  the  cross-wire  system  and  infinite  outside,  which  is  shown  by  the  shaded 
regions.  The  shaded  (blue)  regions  represent  that  these  two  dimensional  region  of  space  is  not 
accessible  for  an  electron.  Thus,  the  motion  of  an  electron  is  confined  within  the  cross  region, 
the  channels  of  which  extend  to  infinity.  The  intersection  region  of  the  cross-wire  system 
is  in  general  a  rhombus  of  side  2b  and  the  width  of  each  nanowires  is  2a.  It  is  well  known 
that  when  the  intersection  geometry  of  the  cross-wire  system  is  a  square,  such  a  classically 
unbound  system  possess  a  quantum  bound  state  [1].  The  motion  of  a  quantum  particle, 


FIG.  22.  Nanowires  in  a  skewed  geometry.  The  potential  is  zero  inside  the  cross-wire  system  and 
infinite  outside.  The  shaded  (blue)  regions  represent  the  regions  not  accessible  for  an  electron.  The 
motion  of  an  electron  is  confined  within  the  cross.  The  intersection  region  of  the  cross- wire  system 
is  a  rhombus  of  side  26,  the  acute  angle  between  its  sides  is  u  and  the  width  of  each  nanowires  is 


2a. 
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such  as  an  electron,  obeys  the  Schrodinger  equation,  as  determined  by  its  Hamiltonian. 
Earlier  literature  sources  (see  Chapter  I),  show  that  such  an  electron  may  either  remain 
localized  in  a  small  region  of  space  (a  bound  state)  or  diffuse  through  the  entire  system  (a 
scattering  state).  We  are  particularly  interested  in  two-dimensional  crossed  nanowires  with 
the  condition  that  they  intersect  at  an  arbitrary  angle.  Previous  studies  on  the  cross-wires 
systems  have  revealed  that  an  electron  in  such  a  system  may  either  be  confined  within  the 
intersection  region  in  a  localized  state  or  decay  along  the  channels  of  the  crossed  wires. 
These  quantum  bound  states  are  responsible  for  the  observed  anomaly  in  the  transmission 
and  reflection  properties  of  crossed  nanowires  systems.  Various  methods  of  observing  these 
states  by  means  of  measurement  of  conductance  were  also  proposed  theoretically  using  two 
dimensional  electron  gas  models.  Furthermore,  much  work  has  been  done  both  in  two  and 
three  dimensional  realistic  system,  finding  new  approach  to  observe  the  effects  of  the  localized 
states  on  the  optical  properties  of  the  material.  We  extend  this  knowledge  by  analyzing  the 
ground  states  binding  energy  of  an  electron  in  the  skewed  geometry,  which  until  now  have 
not  been  studied.  Applying  variational  principle,  we  investigate  conditions  for  a  general 
criterion  for  the  existence  of  a  bound  state  and  the  role  of  tilted  geometry  in  the  bound  state 
energy  of  an  electron.  Also,  the  theory  developed  in  this  work  might  be  useful  to  interpret 
electron  transport  peculiarities  in  realistic  systems  such  as  semiconductor  nanowire  films  and 
carbon  nanotube  bundles.  We  assume  that  our  results  may  have  practical  applications  in 
the  fields  of  materials  science  and  in  the  design  of  quantum  information  devices.  Therefore, 
the  understanding  of  electron  tunneling  properties  in  these  cross-wires  system  is  of  extreme 
relevance  for  device  applications.  The  main  objectives  of  this  thesis  work  can  be  enumerated 
as  follows: 

1.  To  formulate  the  problem  of  cross-wire  system  in  the  skewed  geometry  and  to  study 
angular  dependence  of  the  bound  state  energy  for  an  electron  trapped  at  the  intersec¬ 
tion. 

2.  Use  the  results  of  this  model  to  interpret  the  electron  transport  properties  in  a  realistic 
systems. 

In  the  following  section,  we  describe  the  procedure  to  study  the  cross-wire  system  in  the 
skewed  geometry. 
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B.  Geometry  of  Nanowires 

We  consider  the  coordinate  system  (x,  y)  rigidly  fixed  at  the  origin  (0,  0)  which  is  the 
center  of  the  crossed  nanowires  intersection  (see  figure  23(a)).  We  define  two  coordinate 
systems  attached  with  each  nanowires.  The  coordinate  system  (x\,y{)  describes  the  hor¬ 
izontal  channel  such  that  aq-axis  is  along  the  channel  and  yi-axis  is  perpendicular  to  the 
channel  as  shown  in  the  figure  23(b).  Similarly,  the  coordinate  system  (x2? y2)  is  such  that 
x2-axis  is  along  the  vertical  channel  .  The  origin  of  all  the  coordinate  systems  coincides  but 
are  different  in  their  orientations  and  preserve  right-handedness.  These  coordinate  systems 
are  related  to  the  original  coordinate  system  by  the  following  transformation  relation 


(b)  Vertical  nanowire 


FIG.  23.  The  diagram  (figure  is  highly  exaggerated  for  clarity)  show  the  coordinate  systems  used 
for  the  X-shaped  quantum  nanowires.  The  two  nanowires  (  of  equal  width  2a  )are  shown  separately: 
(a)  The  coordinate  system  rigidly  fixed  to  the  horizontal  nanowire  (channel  along  x-axis)  is  rotated 
counterclockwise  by  an  angle  0\.  (b)  The  coordinate  system  rigidly  fixed  to  the  vertical  nanowire 
(channel  along  y-axis)  is  rotated  counterclockwise  by  an  angle  #2. 


(8) 
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The  above  transformation  relation  describes  the  orientation  of  the  new  axes  attached  to  each 
nanowires  relative  to  the  axis  fixed  at  the  center  of  the  cross.  All  of  these  three  coordinate 
systems  are  orthogonal.  The  unit  vectors  corresponding  to  direction  of  the  axes  of  these 
coordinate  systems  satisfy 

x.y  =  xi.yj  =  x2.y2  =  1  (10) 

Up  to  now  nothing  has  been  done  to  the  geometry  of  the  crossed  nanowires.  The  config¬ 
uration  is  still  the  same.  We  describe  below  a  scheme  to  get  the  geometry  appropriate  to 
our  problem.  To  begin  with,  the  horizontal  channel  is  rotated  in  counter-clockwise  direction 
by  an  angle  6\  with  respect  to  the  original  coordinate  system  (xi,yi).  The  coordinate 
transformation  is  given  by 


(  X>1 1 

cos  9i 

—  sin  9i 

\y,i ) 

sin  9\ 

COS  01 

cos  9\  —  sin  91 
sin  9i  cos  9i 

In  the  next  step,  we  rotate  the  vertical  channel  without  affecting  the  horizontal  one.  Then 
the  coordinate  transformation  for  the  rotation  of  the  vertical  channel  in  counter-clockwise 
direction  by  an  angle  02  relative  to  the  original  coordinate  system  (x2,  y2)  is  given  by 


y  —  cos  9 

It  is  clear  from  the  above  transformation  that 

nr^  —  nr  nr\c  f).  —  n  /  cirt 


sin  02  cos  02  J 

cos  92  —  sin  9-2 
sin  02  cos  02 

sin  92  cos  92 
-  cos  92  sin  92 


x[  =  x  cos  9 1  —  y  sin  9\ ,  y\  —  x  sin 


x'2  =  x  sin  92 


92  +  ycos92,  y'2  = -x  cos 


(12) 


(13) 

(14) 
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Now,  we  have  two  choices  to  select  the  new  coordinate  system.  For  the  present  configuration, 
we  take  {x\ ,  x'2)  as  the  new  basis  with  the  axes  along  x\  and  x'2  direction  which  we  relabel 
them  by  x'  and  y'  respectively.  Hence,  the  nanowires  intersection  geometry  referred  to  the 


FIG.  24.  The  figure  shows  the  intersection  geometry  of  two  nanowires  after  two  successive  or¬ 
thogonal  transformations.  The  intersection  region  is  a  rhombus  of  side  2b.  The  width  (2a)  of  the 
nanowires  remains  constant  in  the  new  geometry. 

new  coordinate  system  (x'^x^)  is  a  rhombus  (see  figure  24).  The  corners  of  this  rhombus 
are  related  to  the  width  of  the  channel  and  the  angle  between  the  axes.  The  new  and 
the  final  coordinate  system  ( x',y' )  is  related  to  the  original  system  through  the  following 
transformation  formula 


x2 


(15) 


xi  \  cos  0 1  —  sin 

sin  62  cos  02 

The  determinant  of  the  transformation  matrix  from  the  unprimed  to  the  primed  set  is  given 
by 


A 


cos  0\  —  sin  0\ 
sin  02  cos  0o 
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=  cos  0 1  cos  02  +  sin  6 1  sin  02 

=  cos(0\  —  62)  =  cos  6  (16) 


The  angular  difference  6  =  (0\  —  62)  gives  the  angle  by  which  one  nanowire  is  rotated 
with  respect  to  the  other.  In  other  words,  this  is  the  measure  of  asymmetry  of  the  crossed 
nanowires  as  shown  in  the  figure  24.  If  two  wires  are  rotated  such  that  0\  =  62  then  they 
intersect  symmetrically  at  right  angle.  If  however  the  condition  0i  —  02  =  n/ 2  is  satisfied, 
then  the  two  wires  completely  overlap  and  form  a  single  wire.  The  inverse  transformation 
from  the  primed  system  to  the  unprimed  system  is  given  by  the  following  transformation 
relation: 


x 

V 


1 

cos  02  sin  0\ 

AT 

A 

—  sin  02  cos  0\ 

UJ 

(17) 


The  Jacobian  of  the  transformation  from  the  unprimed  to  the  primed  set  is  given  by 


J 


d(x,y) 

d(x',y') 


dx 

dx 

dx' 

dy ' 

dy 

dy 

dx1 

dy ' 

cos  02  /A  sin0  i/A 
—  sin02/A  cos^i/A 


1 

A2 

1 

A2 

1 

A2 


cos  02  sin  0i 
—  sin  02  cos  0\ 

(cos  0\  cos  02  +  sin  0\  sin  02) 

cos(01  -  02)  =  ^  (A)  =  ^ 


sec  0 


(18) 


The  first  and  second  partial  derivatives  transform  from  the  unprimed  coordinates  system  to 
the  primed  set  as 
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The  2D  Laplacian  operator  in  the  new  coordinates  system  becomes 


V2  = 


8 2  d 2 

+ 


2  sin(6>1  —  d2 


d2 


=  V72 


dx'2  dy'2  dx'dy' 

The  Hamiltonian  operator  in  the  new  coordinate  system  becomes 


2  sin  6 


d 2 


dx'dy' 


(19) 


H  =  -2L  fv'2  -  2  sin  8-  ° 


(20) 


2m  V  dx'dy'  / 

We  consider  that  the  electron  is  constrained  by  the  potential  which  is  zero  inside  the  chan¬ 
nels  and  infinite  outside  which  is  given  by 


V(x',y')  =  0  for  \x\  <  b  or  \y'\  <  b 
=  oo  elsewhere 

Our  task  is  to  solve  the  free  particle  Schrodinger  equation  for  an  electron  subject  to  the 
above  potential. 

(v'2-2sm^  +  ^)^'!'')  =  0  (21) 
This  equation  clearly  describe  the  behavior  of  an  electron  in  a  two  dimensional  asymmetric 
crossed  nanowires  in  new  geometry.  We  take  Q\  >  02  so  that  sin  9  is  positive.  Here,  the 
width  of  the  nanowire  channels  remain  constant.  This  is  the  unique  problem  where  there 
is  change  in  the  angle  of  the  intersection  of  the  cross-wire  system  without  any  variation 
in  the  width.  Paranjape  calculated  the  ground-state  energy  of  an  electron  trapped  at  the 
intersection  of  a  cross  formed  by  two  quantum  wires  where  the  widths  of  the  wires  forming 
the  cross  were  assumed  to  vary  independently  [13].  He  studied  the  variation  in  the  energy 
with  the  ratio  of  the  widths  and  observed  that  the  energy  decreases  as  one  of  the  widths  is 
allowed  to  increase.  But  our  model  is  based  on  the  angular  variation  rather  than  the  size  of 
the  channels.  In  what  follows,  we  discuss  the  approach  used  to  formulate  the  problem. 


C.  Method  and  Assumption 


1.  Method 

The  approach  used  here  to  study  the  oblique  crossed  wires  system  is  based  on  the  varia¬ 
tional  principle.  This  is  also  known  as  Rayleigh-Ritz  variational  method  which  gives  a  more 
accurate  estimate  of  the  ground  state  energy  than  the  perturbation  theory.  It  is  the  non- 
perturbative  method  to  the  approximate  calculations  of  eigenvalues  [14].  In  what  follows, 
we  give  a  sketch  of  this  method  itself  [15].  The  Schrodinger  equation  can  be  derived  from 
the  variational  principle,  which  requires  that 

j  =  {<f>\H\<i>)  (22) 

is  stationary  with  respect  to  the  independent  variation  of  0  and  0*  subject  to  the  normal¬ 
ization  condition  . 


(m  =  i  (23) 

Using  the  method  of  Lagrange  multiplier,  we  construct  a  function  given  by 

J(e)  =  {<j>\H\<l>)-e{<i>\<j>)  (24) 

where  e  is  the  undetermined  multiplier.  The  requirement  for  the  function  J{e)  to  be  sta¬ 
tionary  leads  to 


5J(e)  =  0 

(5<f>\H  -  e\<f>)  +  (<f>\H  -  e\S(f>)  =  0  (25) 

One  should  be  cautious  in  putting  each  term  in  the  above  equation  equal  to  zero  because 
the  variations  (50|  and  1 8(f))  are  not  independent.  In  general,  we  consider 

S(p  =  5u  +  iSv  (26) 

where  Su  and  Sv  are  arbitrary  real  and  independent  variations.  Then  we  get, 

5J(e)  =  ( 8u\H  —  e\(j))  +  —  e\ 8u)  +  (iSv\H  —  e\4>)  +  (4>\ H  —  e\ iSv) 

=  ( 8u\H  —  e\(f>  +  (f>*)  —  i(5v\H  —  e\ <fi  —  0*) 


(27) 
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For  arbitrary  Su  and  Sv,  one  can  write 

(n-e)  10  +  0*)  =  0 

(h  -  e)  10  -  0*)  =  o  (28) 

On  adding  and  subtracting  these  equation,  one  get 

(£-e)  10)  =  0 

(^-e)|0*)  =  (0|(^-e)=O  (29) 

This  is  equivalent  to  equating  each  term  in  earlier  equation  to  zero.  Thus,  £  has  the  signifi¬ 
cance  of  real  eigenvalue  of  H,  which  is  given  by 

H\<t>)=e\4>),  H\4?)  =  e\4?) 

e  =  =  Wm  =  E  (30) 

Thus,  the  variational  principle  leads  to  a  Schrodinger  equation.  This  means  that  the  solution 
of  variational  problem  is  any  solution  to  the  time  independent  Schrodinger  equation.  In  other 
words,  if  any  variation  is  allowed 

5<0|i0|0)  =  0  (31) 

is  equivalent  to 

#  10)  =  m  (32) 

Suppose  E  is  defined  as  above.  Let  Em  and  (pm  be  an  eigenvalue  and  the  corresponding 
eigenstate  of  H,  respectively.  Then 

5E  =  E-Em 
=  <0|if  |0)  -  Em 
=  <0|iL|0>  -  <0| <j>)Em 

=  (<f>\H  -  Em\(f>)  (33) 

If  0  and  0m  differ  by  50,  that  is 


0  —  0m  +  50 


(34) 
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Then  the  variation  in  energy  SE  becomes 

SE  =  ( 5(j)\H  -  Em\5<f>)  +  O  ((<50)2)  (35) 

It  means  that  to  the  first  order  of  approximation, 

SE  =  0  (36) 

This  is  variational  principle.  Thus,  the  best  approximation  to  Ern  is  obtained  by  varying  0, 
or  the  parameters  on  which  it  depends  so  that  above  relation  is  obeyed.  The  Hamiltonian 
operator  H  which  is  hcrmitian,  generates  a  complete  orthonormal  set  of  eigenfunctions: 

H\4>m)  =  E\4>m),  (0m|0n)  =  Smn  (37) 


We  now  expand  the  the  trial  function  <f>  in  terms  of  complete  orthonormal  set  of  eigenfunc¬ 
tions 


£' 


=  >  a. 


The  energy  or  the  expectation  value  of  the  Hamiltonian  is  given  by 

E  \(f>)  ^  ^ (Qm0m | H\ 

m,n 

=  ^  ^  OjmClnEn5mn  =  ^  ^  |Om|  Em 

m,n  m 

and 


(38) 


(39) 


1  (010)  ^  ^  1 0n)  ^  y  UmUn(5mn  ^  ^  (lm  (40) 

m,n  m,n  m 

If  Eq  is  the  lowest  energy  among  the  spectrum  Em,  then  we  have 

E  >  E0  |  cim  | 2 

m 

Therefore,  E  >  E0  (41) 

In  order  to  obtain  a  upper  bound  for  some  kth  level  Ek  of  the  spectrum,  we  select  a  <f> 

which  is  orthogonal  to  the  eigenfunctions  0O,  0i, . ,4>k- 1  for  all  the  lower  levels.  Then  the 

coefficients  cii,a2 . ,  l  are  all  zero.  The  lowest  energy  which  now  occurs  for  E  is  Ek. 

It  follows  that  E  >  Ek. 

In  atomic  and  molecular  physics,  most  methods  for  finding  the  energy  levels  are  varia¬ 
tional  schemes  [16].  The  Hartree  method  of  self-consistent  fields  which  is  one  of  the  ways 
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to  find  the  central  field  is  the  generalization  of  the  variational  principle.  In  many-electron 
atoms,  the  trial  function  (j)  is  expressed  as  a  linear  combination  of  known  functions  ut  as 

r 

0  =  Y  °iUi  (42) 

i=  1 

where  c\s  are  the  variational  parameters.  The  functions  tq  are  not  necessarily  orthogonal  so 
that  we  have  the  non-zero  overlap  integral  given  by 


Sij  \  U-i  I  U  )  ) 


(43) 


The  expectation  value  of  energy  E  for  the  above  choice  of  the  trial  function  is  given  by 

(0i  m 


E  = 

E  = 


(010) 

^2ij  Ci '  ./  H'  j 


^2ij  Ci  CjSij 

E  CiCjSij  =  CiCjHij 
ij  ij 

where  //,?  =  (uj\H\uj). 

Differentiating  both  sides  with  respect  to  one  of  the  coefficients,  say  c*k,  we  get 

Ci  CjSij  +  E  CjSkj  =  CjHkj 
k  ij  j  j 

The  minimum  energy  E  is  obtained  by  putting 

d  E 


(44) 


(45) 


del 


=  0 


or 

r 

Y.  Cj  ( Hkj  -  ESkj)  =  0,  k  =  1,  2 . ,  r  (46) 

1=1 

The  corresponding  secular  equation  is 

=  0  (47) 

Hence,  with  the  given  choice  of  the  trial  functions  the  matrix  elements  Hij  of  the  Hamiltonian 
and  the  overlap  integrals  St]  are  evaluated.  If  an  orthonormal  basis  is  used,  the  above  secular 


Hu  -ES 


li 


Hr  1  -  ES, 


r  1 


H\r  —  ES 


1  r 


Hr 


ESr 
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equation  is  greatly  simplified  because  =  Sij.  In  this  case,  the  secular  determinant  becomes 


H 


11 


E 


Hrl  —  E 


H 


1  r 


Hr. 


E 


E 


=  0 


(48) 


In  either  case,  the  secular  determinant  for  r  basis  functions  give  rth  order  polynomial  in  E 
which  is  solved  for  r  different  roots,  each  of  which  approximates  a  different  eigenvalue.  The 
r  roots  are  the  possible  values  of  E,  the  lowest  one  being  an  upper  bound  to  the  ground 
state  energy  of  the  molecular  system.  The  variational  method  is  a  basis  for  Hartree-Fock 
theory  and  the  configuration  interaction  method  used  in  the  electronic  structure  of  atoms 
and  molecules. 

In  the  variational  method,  a  fixed  H  is  taken  and  it  can  be  shown  how  the  matrix 
element  of  H  varies  as  the  trial  function  is  modified.  There  is  a  theorem  after  Hellmann  and 
Feynmann  which  shows  how  the  energy  of  a  system  varies  when  the  Hamiltonian  if,  which 
depends  on  a  parameter  f3,  changes,  i.e., 


dE  _  dH 
~d/3  ~ 


(49) 


The  significance  of  this  theorem  is  that  the  operator  ^  may  be  simple  in  its  form.  As 
an  example,  if  H  —  H0  +  /3x,  then  ^  =  x,  which  does  not  involve  H0,  the  unperturbed 
Hamiltonian.  Hence,  using  trusted  trial  functions  can  be  evaluated  from  which  we  obtain 
the  energy  eigenvalues  of  the  system. 

In  practice,  if  we  have  a  trial  function  (j)  as  a  function  of  some  variational  parameter  a, 
then  the  expectation  value  of  the  Hamiltonian  is  given  by 


,fr(  nv  =  {(p(r,a)\H\(j)(r,a)) 
[  [  )}  (0(r,a)|0(r ,«)) 


(50) 


The  ground  state  energy  is  equal  to  the  value  of  this  ratio  evaluated  using  variational  param¬ 
eter  (a)  that  minimizes  this  ratio.  It  just  involves  the  evaluation  of  two  integrals  without 
having  to  solve  any  differential  equation.  The  variation  theorem  states  that  E  >  Eq  for  any 
trial  function  (f>.  The  equality  holds  if  the  trial  function  (j)  is  the  true  ground  state  wave 
function.  A  trial  function  is  set  up  in  a  flexible  form  in  terms  of  parameters  making  a  guess. 
Its  parameters  are  varied  to  reduce  the  trial  energy  E  to  its  lowest  value  possible.  The  result 
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is  an  upper  bound  to  the  ground  state  energy.  The  energy  approaches  the  true  ground  state 
energy  from  the  top.  This  means  that  the  energy  obtained  from  the  variational  method  is 
always  above  the  true  ground  state  energy. 


2.  Assumption 


We  assume  that  the  shape  of  the  potential  which  confines  the  electron  to  the  channel 
is  not  significant  in  determining  the  existence  of  bound  state  at  the  intersection  of  the  two 
wires  [1].  This  is  because  only  the  corners  are  responsible  for  the  existence  of  bound  state 
[9].  It  has  been  proved  that  in  a  single  straight  channel  bound  state  does  not  exist  due 
to  absence  of  corners  and  bents.  Thus,  the  potential  has  nothing  to  do  with  the  electron 
confinement.  In  general,  potential  can  be  taken  to  be  constant,  harmonic  oscillator  type 
(parabolic)  or  any  complex  function.  Hence,  it  is  sufficient  to  take  the  potential  to  be  zero 
within  the  rhombus  region  and  the  semi-infinite  channels.  The  region  outside  the  cross-wire 
system  is  inaccessible  to  an  electron  as  the  walls  does  not  allow  tunneling.  The  walls  are 
rigid  and  impenetrable.  We  choose  proper  trial  functions  and  use  national  principle  to  find 
the  upper  bound  to  the  ground  state  energy  of  the  system.  The  electron  is  constrained  in  a 
X-shaped  channels  in  the  primed  coordinate,  so  the  trial  wave  functions  are  taken  as 


(l  -  e-«  for 

^1  —  ^  j  e  ^  for 
^1  —  e~~^  for 


\x'\  <  b,  \y'\  <  b 
\x'\  >  b,  \y'\  <  b 
\y'\  >  b ,  \x'\  <  b 


0 


otherwise 


(51) 


where  A  is  the  normalization  constant  which  may  depend  on  the  variational  parameter  (a) 
and  also  on  the  asymmetry  of  the  wires.  These  trial  functions  are  such  that  they  vanish 
outside  the  rhombus  of  side  2b.  The  origin  of  the  coordinate  system  y')  is  taken  at  the 
center  of  the  rhombus.  The  region  of  localization  is  a  rhombus  with  the  corners  at  (±6, ±6), 
where  b  A  =  a  and  w  =  2a  is  the  constant  width  of  the  nanowires  as  shown  in  figure  24.  For 
simplicity,  we  write  the  relative  rotation  and  the  Jacobian  respectively  as 

9\  —  02  =  9  and  A-1  =  sec  9 


In  the  present  work  that  follows,  the  above  relations  have  been  used  extensively  to  avoid 
computational  complexity. 
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III.  TRIAL  FUNCTIONS 


A.  Trial  functions  for  the  skewed  cross-wire  geometry 


The  shape  of  the  potential  is  insignificant  in  determining  the  bound  states  in  the  cross¬ 
wire  system.  The  potential  is  zero  within  the  two  dimensional  channel  and  infinite  outside 
it.  The  electron  is  located  inside  the  crossed  intersection  with  both  channels  extending  to 
infinity.  An  electron  is  in  the  region  of  zero  potential  so  that  no  external  force  act  on  it.  It 
has  total  energy  only  dne  to  its  two  components  of  momentum  (horizontal  and  the  transverse 
component).  It  will  be  clear  later  in  the  chapter  that  when  an  electron  is  inside  the  channel 
far  away  from  the  cross  region,  the  corresponding  transverse  component  of  the  momentum 
is  quantized.  This  is  a  typical  property  of  a  quantum  wire.  This  feature  drastically  changes 
within  the  region  of  intersection  of  two  identical  wires.  The  crossed  intersection  behaves  as 
two  dimensional  potential  well.  The  electron  can  be  trapped  and  bound  to  the  crossed  region 
without  any  further  decay  into  the  continuum  of  the  channel.  The  electron  is  locally  bound 
in  the  crossed  wires  which  is  inherently  quantum  in  origin.  The  trapped  state  manifest 
itself  as  quantum  dot  in  two  dimension.  The  size  of  quantum  dot  is  determined  by  the  area 
of  intersection  and  the  nature  of  the  corners  (sharp  edge  or  rounded).  The  trial  functions 
is  selected  making  a  guess  and  satisfying  some  conditions.  The  trial  functions  within  the 
intersection  region  form  a  roof  with  rooflines  and  falls  exponentially  just  outside  the  rhombus. 


<$>(x',y')=A{ 


\x'\  <  b,  \y'\  <  b 
\x'\  >  b,  \y'\  <  b 
\y'\  >  b,  \x'\  <  b 


0 


otherwise 


(52) 


where  A  is  the  normalization  constant  which  may  depend  on  the  variational  parameter  (a) 
and  also  on  the  asymmetry  of  the  wires.  These  trial  functions  are  such  that  they  vanish 
outside  the  rhombus  of  side  2b.  The  coordinate  system  ( x',y' )  is  taken  at  the  center  of  the 
rhombus.  The  region  of  localization  is  a  rhombus  with  the  corners  at  (±6, ±6),  where  bA  =  a 
and  w  =  2 a  is  the  constant  width  of  the  nanowires  as  shown  in  the  figure  24.  Using  the 
coordinate  transformation  from  Chapter  2,  the  three  dimensional  plots  for  typical  values  of 
rotation  angles  (0)  such  as  0°,  30°,  45°,  60°  and  90°  are  neatly  displayed.  These  plots  show 
the  behavior  of  trial  functions  used  in  the  calculations. 
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FIG.  25.  Trial  function  for  crossed  wire  system  without  any  rotation  (9  =  0), 
at  right  angle. 


the  wires  intersect 


FIG.  26.  Probability  density  for  crossed  wire  system  without  any  rotation  ( 9  =  0),  the  wires 
intersect  at  right  angle. 


Figure  25  shows  the  trial  function  without  any  rotation  of  the  crossed  wire  system.  These 
wires  intersect  at  right  angle  and  the  region  of  intersection  is  a  square  geometry.  The 
rooflines  cross  each  other  at  right  angle.  This  configuration  has  the  trial  functions  decaying 
exponentially  along  the  channels  of  the  cross  and  vanishing  abruptly  at  the  corners.  This 
feature  is  more  clear  from  the  probability  density  plot  [Fig.26].  The  spatial  extend  of  the 
trial  wave  function  is  localized  well  inside  the  intersection  region  and  also  the  probability 
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density  is  sharply  peaked.  This  feature  make  the  system  more  likely  to  be  stable  and  in  fact 
it  is  the  most  stable  state. 


FIG.  27.  Trial  function  for  rotation  of  the  axes  by  angle  9  = 


vr/6. 


FIG.  28.  Probability  density  for  rotation  of  the  axes  by  angle  8  =  7r/6. 

Figure  [Fig. 27]  shows  the  plot  of  trial  function  for  the  crossed  wires  in  which  one  of  the 
wires  is  rotated  by  30°.  These  wires  intersect  at  an  angle  of  60°  and  the  region  of  intersection 
is  a  rhombus.  The  rooflines  cross  each  other  at  an  angle  60°.  This  configuration  has  the 
trial  functions  decaying  exponentially  along  the  channels  of  the  cross  and  vanishing  at  the 
corners.  This  feature  is  more  clear  from  the  probability  density  plot  [Fig.28].  The  spatial 
extend  of  the  trial  wave  function  is  little  bit  broaden  outside  of  the  intersection  region  and 
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also  the  probability  density  is  less  sharply  peaked.  The  system  is  still  stable  and  not  like 
that  of  the  unperturbed  geometry. 


FIG.  29.  Trial  function  for  rotation  of  the  axes  by  angle  6  =  n/4. 


FIG.  30.  Probability  density  for  rotation  of  the  axes  by  angle  9  =  7t/4. 

Figure  [Fig. 29]  shows  the  plot  of  trial  function  for  the  crossed  wires  in  which  one  of  the 
wires  is  rotated  by  45°.  These  wires  intersect  at  an  angle  of  45°  and  the  region  of  intersection 
is  a  rhombus.  The  rooflines  cross  each  other  at  an  angle  45°.  This  configuration  has  the 
trial  functions  decaying  exponentially  along  the  channels  of  the  cross  and  vanishing  at  the 
corners.  This  feature  is  more  clear  from  the  probability  density  plot  [Fig.30].  The  spread  of 
the  wave  function  outside  the  intersection  region  is  more  due  to  electron  tunneling  and  is 
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less  stable  but  still  it  is  a  bound  state. 


FIG.  31.  Trial  function  for  rotation  of  the  axes  by  angle  6  =  7t/3. 


FIG.  32.  Probability  density  for  rotation  of  the  axes  by  angle  8  =  7r/3. 

Figure  [Fig. 31]  shows  the  plot  of  trial  function  for  the  crossed  wires  in  which  one  of  the 
wires  is  rotated  by  60°.  These  wires  intersect  at  an  angle  of  30°  and  the  region  of  intersection 
is  a  rhombus.  The  rooflines  cross  each  other  at  an  angle  30°.  This  configuration  has  the 
trial  functions  decaying  exponentially  along  the  channels  of  the  cross  and  vanishing  at  the 
corners.  This  feature  is  more  clear  from  the  probability  density  plot  [Fig.32],  The  trial  wave 
function  is  spatially  broaden  due  to  electron  tunneling  and  the  stability  of  the  system  is 
lowered. 
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FIG.  33.  Trial  function  for  rotation  of  the  axes  by  angle  0  =  7t/2. 


FIG.  34.  Probability  density  for  rotation  of  the  axes  by  angle  9  =  tt/2. 

Figure  [Fig. 33]  shows  the  trial  function  for  the  crossed  wires  in  which  one  of  the  wires 
is  rotated  by  90°.  These  wires  merge  into  a  single  channel.  The  two  channels  completely 
overlap  each  other.  There  are  no  rooflines.  This  feature  is  more  clear  from  the  probability 
density  plot  [Fig. 34].  It  will  be  shown  later  that  this  geometry  is  not  favorable  to  form  a 
quantum  bound  state.  The  electron  wave  function  completely  spread  out  along  the  channels. 
The  probability  density  show  that  an  electron  can  exist  anywhere  along  the  channel.  This 
is  the  case  of  free  propagation  of  an  electron  along  the  channels  (scattering  states). 
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B.  Localization  regions  of  an  electron 

The  electron  can  be  located  inside  the  channel  extending  to  infinity  in  both  direction  along 
x'  and  y'  axes.  The  geometry  of  the  nanowires  intersection  is  shown  in  the  figure  35.  This 
figure  is  highly  exaggerated  so  that  the  different  regions  can  be  visualized.  Also,  evaluation  of 
the  normalization  constant  and  the  matrix  elements  of  the  Hamiltonian  requires  the  regions 
to  be  properly  specified.  It  is  for  this  reason,  the  regions  are  partitioned  accordingly.  The 
intersection  region  is  in  general  a  rhombus  with  the  corners  at  (±6,  ±6).  The  side  of  this 
rhombus  is  2b.  The  width  (w  =  2a)  of  the  nanowire  is  same  for  both  the  channels.  This 
is  because  orthogonal  transformation  to  the  right-handed  coordinate  system  was  used  to 
obtain  this  oblique  geometry  which  is  our  model  system  as  discussed  in  the  objective  of 
this  work.  And  under  such  transformation,  width  remains  invariant.  The  spatial  extend 


FIG.  35.  Localization  regions  of  an  electron  in  the  oblique  cross-wire  system.  6  is  the  angle  of 
rotation  of  the  nanowire  and  a  is  the  semi-width  of  the  nanowires.  The  corners  of  the  rhombus  are 
(±6,  ±6)  referred  to  the  primed  coordinate  system. 
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of  the  wave  functions  show  that  the  intersection  has  the  greatest  probability  of  electron 
localization  and  the  probability  decreases  as  we  go  along  the  arms  of  the  cross-wire  system. 
In  general,  there  are  two  regions  of  localization  within  the  cross-system.  The  electron  may 
be  localized  within  the  intersection  region  (rhombus)  and  in  the  channels.  The  rhombus 
can  be  divided  into  four  parts  labeled  as  R\,R2,Rg  and  R 4  located  in  the  first,  second, 
third  and  fourth  oblique  quadrants  respectively.  There  are  four  semi-infinite  channels  in  the 
X-shaped  geometry.  Each  of  the  four  channels  (arms)  extending  to  infinity  are  divided  into 
two  portions  making  eight  channels.  Hence,  there  are  eight  semi-infinite  strips  labeled  as 
Ci,  These  regions  have  different  geometrical  orientation  with  respect  to  the  origin  of 

the  oblique  coordinate  system.  This  means  that  the  boundaries  in  the  two  dimensional  space 
are  different  for  these  regions.  The  angle  6  represent  the  relative  rotation  of  the  channels 
of  the  cross-wire  system.  This  fact  is  clear  from  the  geometry  of  the  nanowires  discussed 
in  Chapter  11(B).  This  rotation  angle  is  different  from  the  angle  of  intersection  (a;)  of  the 
cross-wire  system  (see  figure  22).  In  the  present  work,  all  the  calculations  are  carried  out 
in  terms  of  the  rotation  angle  6  rather  than  the  angle  of  the  intersection  of  the  cross-wire 
system  (a;). 


C.  Normalization  of  the  trial  functions 


To  obtain  the  normalization  constant  A  appearing  in  the  sets  of  trial  functions,  the 
normalization  condition  in  the  unprimed  coordinate  system  can  be  written  as 


($(x,y)\$(x,y)) 


&*(x,y)$(x,y)dxdy  =  1 


Sl&(x,y) 


(53) 


In  the  new  coordinate  system,  this  condition  becomes 


Qe(x'  ,y') 


where,  is  the  domain  of  integration. 

In  practice,  it  is  not  required  to  evaluate  the  Jacobian  for  the  orthogonal  transformation 
of  coordinate  system.  The  Jacobian  of  the  coordinate  transformation  is  a  constant  (see 
Eqn.18).  Since,  variational  principle  is  used,  it  will  appear  in  both  the  numerator  and 
denominator  of  the  Rayleigh-Ritz  ratio  (see  Eqn.50),  so  that  it  can  be  ignored.  Hence,  it  is 
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independent  of  the  coordinate  system  selected  for  the  calculations  as  long  as  the  Jacobian 
is  constant.  But  in  the  present  context,  Jacobian  is  included  for  convenience  even  though  it 
does  not  alter  the  final  result.  The  integral  is  evaluated  over  X-shaped  planar  regions  which 
represent  the  two  dimensional  spatial  extend  for  an  electron  described  by  the  sets  of  trials 
wave  functions  discussed  earlier. 

In  order  to  find  the  normalization  constant,  we  divide  the  region  into  four  oblique  quad¬ 
rants  represented  by  Ri,R2,  R3  and  and  eight  channels  labeled  as  C\ ,  C2 ,  C3 ,  C4 ,  C5 ,  C6 ,  C7 
and  Cg  as  shown  in  the  figure  24.  The  symmetry  of  the  trial  functions  in  different  regions  of 
the  nanowire  intersection  can  be  used  to  evaluate  these  integrals  explicitly.  In  what  follows, 
the  steps  taken  to  arrive  at  the  normalization  constant  will  be  described  for  the  given  sets 
of  trial  wave  functions.  For  detailed  calculations  of  the  integrals  refer  to  the  Appendix  I.  In 
order  to  increase  the  efficiency  of  computation,  some  integrals  were  evaluated  using  math¬ 
ematical  handbook  [17,  18].  The  values  of  the  integrals  for  the  oblique  quadrants  defined 
by  Ri,R2,  R-.i  and  i?4  are  as  follows: 

R1  =  R2  =  R3  =  Ra  =  11  A2  b2  sec  9  e~2a  (55) 

18 

The  total  contribution  to  the  normalization  integral  from  the  regions  Rn  is  given  by 

4 

h  =  Y,Ri  =  ±Ri 

i=l 

=  4  A2b2  sec  9  e~2a 

=  ^  A2  b2  sec  9  e~2a  (56) 

9 

and  the  values  of  the  integrals  for  eight  channels  labeled  as  C4,  C2,  C3,  C4,  C5,  Cq,  C'7  and  Cg 
are  as  follows: 

Ci  =  C2  =  C3  =  C4  =  C5  =  Ce  =  C7  =  C8=^-b2  sec 9 e~2a  (57) 

6a 

The  total  contribution  to  the  normalization  integral  from  the  channels  Cn  is  given  by 

8 

Ji  =  J2Ci  =  8Ci 

%=i 

=  8  (—  b2  sec  9e~2a\ 

\6a  J 

4  A2 

^  ^  ,2  n  —2a 

=  - b  sec  9  e 

3  a 


(58) 
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Hence,  from  the  normalization  condition,  we  have 


h  +  Ji  — 


22 

~9~ 


A2b 2  sec  9e~2a  +  —b2 
3  a 

2  A2  b2  sec  9  e~2a 
9  a 


sec  9  e~2a  = 
(11a  +  6)  = 


1 

1 

1 


A2 


9ae2a 

2  b2  (11a  +  6)  sec  9 


A 


9  ae2a 

(11a  +  6)  sec  9 


(59) 


This  is  the  required  normalization  constant  for  the  given  set  of  trial  functions.  The  normal¬ 
ized  trial  functions  becomes 


$(x',y')  = 


9a  e 


2a 


2  b 2  (11a  +  6)  sec  9 


i  _  WjA 
1  b2 


i  _  M 

1  b 


e  “  for 

_  OL  \x'  | 


x'\  <  b,  \y'\  <  b 
for  |x7|  >  6,  \y'\  <  b 

^1  —  ^  j  e~ for  \y'\  >  b,  \x'\  <  b 


(60) 


otherwise 


In  the  following  Chapters,  we  retain  the  original  normalization  constant  as  A  in  the  trial 
functions  to  carry  out  the  calculations  in  order  avoid  the  mathematically  complicated  ex¬ 
pression.  The  value  of  the  normalization  constant  evaluated  above  is  used  explicitly  after 
the  final  results  of  the  calculations. 
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IV.  EXPECTATION  VALUE  OF  THE  HAMILTONIAN 


Channel  matrix  elements 


The  skewed  geometry  is  a  cross  with  the  intersection  at  an  arbitrary  angle  enclosing  a 
region  which  is  a  rhombus  of  side  2b.  The  arms  of  this  rhombus  extend  to  infinity  in  both 
horizontal  and  vertical  directions.  These  arms  are  called  channels  of  the  cross-wire  system 
(see  figure  35).  These  channels  are  labeled  as  C\ ,  C2,  C'3,  C'4,  C5j  C6,  C7  and  U8.  Although, 
there  are  only  four  channels  but  for  the  computational  convenience,  each  channel  is  taken  as 
two  portions  making  eight  channels.  Each  of  these  channels  are  identical  in  width  but  are 
different  with  respect  to  the  orientation  relative  to  the  origin.  Also,  calculation  of  matrix 
elements  of  the  Hamiltonian  for  the  channels  require  boundaries  of  the  area  enclosed  by  the 
channels.  These  features  makes  the  channels  completely  different  with  regard  to  computa¬ 
tion.  Hence,  not  all  the  matrix  elements  are  equal.  The  matrix  elements  of  the  Hamiltonian 
corresponding  to  these  eight  arms  excluding  the  rhombus  are  called  channel  matrix  elements. 
For  detailed  calculations  of  the  channel  matrix  elements  refer  to  the  Appendix  II  (A).  Also, 
to  improve  the  efficiency  of  computation,  most  channel  matrix  elements  were  evaluated  us¬ 
ing  mathematical  handbook  [18].  The  channel  matrix  elements  Jcn  corresponding  to  the 
channels  Cn  (see  figure  35)  are  given  by 

— K2A2e~2a  sec9  \< 


JCI  —  JC2  ~  Jc5  ~  JC6  ~ 


4  m 


-  —  sin  0 
3 


(61) 


and 


J Ci  —  -Am  —  Jr 77  —  J('m  — 


-h2A2e 


-2  a 


sec  9 


>C4 


c  7 


>C  8 


4m 


a  ■  n 

— b  sm  6 
3 


(62) 


The  total  channel  matrix  elements  for  the  channels  Cn  is  given  by 

(H)  c  =  Jci  +  Jc2  +  J C'3  +  JcA  +  Jc5  +  Jc6  +  J C7  +  Jc8 
=  2  JCi  +  2Jc3  +  2  Jc^  +  2Jq7 


=  4  (Jci  +  Jcz) 

^  /  —H2A2e~2a  sec  9  \a 


4  m 

—2 h2A2ae~2a  sec  9 
3  m 


- sin  6 

3 


—h2A2e  2  seed  ra  . 

b - -  +  sm  9 

4  m  L  3 


(63) 
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This  is  the  part  of  the  ground  state  energy  of  an  electron  associated  with  the  channels  only. 
In  the  next  section,  we  calculate  the  part  of  the  ground  state  energy  contributed  by  the 
rhombus  (intersection  region). 

B.  Rhombus  matrix  elements 


The  skewed  geometry  of  the  nanowires  intersection  is  a  rhombus.  This  rhombus  consist 
of  four  identical  segments  in  x'  —  y'  plane.  These  four  oblique  quadrants  are  labeled  as 
Ri,  R-2,  R:i  and  R±  (see  figure  35).  The  geometrical  area  occupied  by  these  quadrants  are  the 
same  but  are  different  with  respect  to  the  orientation  relative  to  the  origin.  Hence,  not  all 
the  matrix  elements  of  the  Hamiltonian  are  equal.  The  matrix  element  of  the  Hamiltonian 
operator  evaluated  over  these  oblique  quadrants  excluding  four  semi-infinite  channels  are 
called  rhombus  matrix  elements.  For  detailed  calculations  of  the  rhombus  matrix  elements 
refer  to  the  Appendix  II  (B).  Also,  to  improve  the  efficiency  of  computation,  most  rhombus 
matrix  elements  were  evaluated  using  mathematical  handbook  [17].  The  rhombus  matrix 
elements  Irtn  corresponding  to  the  oblique  quadrants  Rn  (see  figure  35)  are  given  by 

3  H2  A2  .  „  „  n„, 


In,  — 


sin  9  sec  9  e  =  If 


3h2  A2  _2a 

IR2  = - sm#  sec  de  =  In 

4  m 


The  total  rhombus  matrix  elements  is  given  by 


(H)r  —  Iri  +  Ir2  +  I Rz  +  It 


Hence,  the  rhombus  matrix  elements  add  to  zero.  This  means  that  there  is  no  contribution 
to  the  ground  state  energy  for  an  electron  from  the  intersection  region.  It  will  be  seen  shortly 
in  the  next  section,  that  the  major  contribution  to  the  energy  comes  from  the  rooflines. 


C.  Rooflines  matrix  elements 


The  most  important  feature  of  this  problem  is  the  evaluation  of  the  rooflines  matrix 
elements.  The  rooflines  are  the  regions  of  discontinuity  of  the  trial  functions.  Referring  to 
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the  figure  35,  the  rooflines  are  categorized  into  three  basic  types.  The  first  type  are  the 
outer  four  arms  which  are  defined  by  the  following  relations: 


x1  =  b 

to 

x'  —  oo  and  x'  =  — oo 

to 

x'  —  —  b  |  at 

y'  =  0 

(67) 

y'  =  b 

to 

y’  —  oo  and  y'  =  —  oo 

to 

y'  =  -b  |  at 

x'  =  0 

(68) 

The  first  relation  represents  the  outer  right  and  outer  left  arms  respectively.  The  second 
relation  represents  the  outer  top  and  the  outer  bottom  arms  respectively  (see  figure  35).  The 
Hamiltonian  operator  consist  of  the  kinetic  energy  operator  in  one  direction  (either  x'  or  if 
direction)  and  the  second  order  cross  derivative  term  associated  with  the  skewed  geometry 
only.  For  detailed  calculations  of  the  rooflines  matrix  elements  refer  to  the  Appendix  III 
(A).  Also,  mathematical  handbook  [17]  was  extensively  used  to  improve  the  efficiency  of 
computation.  The  matrix  elements  of  the  Hamiltonian  ( Kn,n  =  1,4)  corresponding  to  the 
outer  right  arm  is  given  by 


A'i  = 


K 2  A2  e  2a  sec  9 
2  ma 


(69) 


Following  the  same  steps  of  calculations  for  the  other  three  arms  with  proper  Hamiltonian 
and  the  trial  functions,  we  get 


I<2  =  K3  =  IU  =  Ki 


(70) 


The  total  contribution  to  the  matrix  elements  of  the  Hamiltonian  from  the  four  outer  arms 
is  given  by 


=  Y.K>  =  4K' 

i—  1 

_  2h2  A2e~2asecd 
m  a 


(71) 


The  second  type  of  rooflines  are  the  inner  four  arms,  which  are  defined  by  the  following 
relations: 


x'  =  0 

to 

x'  —  b 

and 

x'  =  —b 

to  x'  —  0  |  at 

y'  =  o 

(72) 

y'  =  o 

to 

y'  =  b 

and 

y'  —  —b 

to  y'  =  0  |  at 

x  =  0 

(73) 

The  first  relation  represents  the  inner  right  and  inner  left  arms  respectively.  The  second 
relation  represents  the  inner  top  and  the  inner  bottom  arms  respectively  (see  figure  35). 
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Here  too,  the  Hamiltonian  operator  consist  of  the  kinetic  energy  operator  in  one  direction 
(either  x'  or  y'  direction)  and  the  second  order  cross  derivative  term  associated  with  the 
skewed  geometry  only.  For  detailed  calculations  of  the  rooflines  matrix  elements  refer  to  the 
Appendix  III  (B).  Most  of  the  integrals  were  evaluated  using  mathematical  handbook  [17]. 
The  matrix  elements  of  the  Hamiltonian  ( Ln ,  n  =  1,4)  corresponding  to  the  inner  right  arm 
is  given  by 


L  i  = 


h2A2e  2a  sec  9 


2m 


(74) 


With  simple  analogy  with  Li  (  matrix  elements  of  the  Hamiltonian  for  the  inner  right  arm), 
the  matrix  elements  for  the  inner  left,  inner  top  and  the  inner  bottom  arms  are  respectively 
given  by 


L2  —  L\ ,  T3  —  L  \  and  A  4  —  L\ 


(75) 


K„  -  V  /..  -  I  /. 

i=  1 

2  h2  A2  e~2a  sec  9 


m 


(76) 


The  third  type  of  rooflines  are  the  outer  eight  sides,  which  are  defined  by  the  following 
relations: 


y'  =  0 

to 

y'  —  b  and 

y'  =  -b 

to 

y'  =  0  j  at 

x'  =  ±b 

(77) 

x'  =  0 

to 

x'  —  b  and 

x'  =  —b 

to 

x'  =  0  |  at 

y'  =  ±b 

(78) 

Also,  the  Hamiltonian  operator  consist  of  the  kinetic  energy  operator  in  one  direction  (either 
x'  or  y'  direction)  and  the  second  order  cross  derivative  term  associated  with  the  skewed 
geometry  only.  For  detailed  calculations  of  the  rooflines  matrix  elements  refer  to  the  Ap¬ 
pendix  III  (C).  Most  of  the  integrals  were  evaluated  using  mathematical  handbook  [17]. 
Also,  Mathematica  package  was  used  for  the  evaluation  of  the  integrals.  The  integrals  in¬ 
volving  9  and  5  functions  were  manually  calculated.  The  matrix  elements  of  the  Hamiltonian 
(Sn,  n  =  1,  8)  corresponding  to  the  outer  first  side  is  given  by 

h2A2e~2asec9  f  2a  —  1' 


Si  = 


(79) 


2  m  \  6 

Use  of  appropriate  trial  functions  for  the  remaining  seven  sides,  it  can  be  shown  that  the 
matrix  elements  Sn  are  equal.  Hence,  the  total  contribution  to  the  matrix  elements  of  the 
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Hamiltonian  from  the  eight  sides  is  given  by 


Km  =  8  Si  = 


h2A2e  2asec  9 

(2a-l\ 

2m 

l  6  )\ 

2h2A2e~2a  sec  9 
3  m 


(2a  -  1) 


(80) 


D.  Matrix  element  of  the  Hamiltonian 


The  contribution  to  the  matrix  element  of  the  Hamiltonian  for  all  the  rooflines  is  given 


by 


(H^Rooflines-Kj  +  Ku  +  Kju 

2h2A2e~2a  sec  9  2h2A2e~2a  sec  9 
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1  _  2a- 
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2 h2A2e~2a  sec  9 
3  m 


(2a  -  1) 
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(81) 


Hence,  the  matrix  element  of  the  Hamiltonian  for  the  entire  cross-wire  system  consists  of 
the  contribution  coming  from  the  channels,  oblique  quadrants  and  the  rooflines. 


<• H(a ))  =  (H)c  +  (H)r  + 
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mb 2  V  11a  +  6 
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(82) 


mb 2  (11a  +  6) 

where,  the  value  of  the  normalization  constant  A  have  been  used  to  obtain  the  final  expres¬ 
sion.  This  is  the  expectation  value  of  the  Hamiltonian  for  an  electron  in  the  entire  cross-wire 
system.  In  the  following  chapter,  this  expression  is  used  to  evaluate  the  ground  state  energy 
of  an  electron  in  the  oblique  cross-wire  geometry. 
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V.  BOUND  STATE  ENERGY 

A.  Threshold  energy  vs  bound  state  energy 

The  threshold  energy  or  the  propagation  threshold  for  the  channels  intersecting  at  an 
arbitrary  angle  can  be  obtained  by  solving  2D  Schrodinger  equation  given  by 

/  . . d2  2  mE\ .  . 


V'2  -  2  8111(0!  -  e2 


dx'dy'  h2 


$(x',y')  =  0 


or'  (s5  +  h  - 2  sm{e 1  -  #2)  aSv  +  AT  )  Hxl’ y,)  =  ° 


with  the  boundary  conditions  that 


$(±6, y')  —  0  for  \y'\  >  b 
&(x',±b)  =  0  for  \x'\  >b 


Using  the  inverse  transformation 


1  cos  9-2  sin  0]  I  x' 
^  —  sin  02  cos  0i  \  y' 


the  2D  Schrodinger  equation  becomes 

(  d2  d2  2  m  E 


3x2  +  3y2+  h2  0 


with  the  following  new  boundary  conditions 

$(±a,2/)  =  0  for  \y\  >  a 
$(x,  ±a)  =  0  for  |x|  >  a. 

Using  the  technique  of  separation  of  variables,  we  can  write 


From  Eqn.86,  we  get 


d >{x,y)=X{x)Y{y) 


1  d2X(x)  1  d2Y(y)  2rnE 
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1  d2Y{y) 
Y(y )  dy2 


=  —  K, 


=  — 


(91) 
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Hence,  energy  E  is  related  to  the  separation  constants  kx  and  Ky  by 

E  =  L  U  +  d)  <92> 

Since  all  channels  are  equivalent,  we  consider  the  channel  in  the  first  quadrant  such  that 
the  wave  is  guided  along  x-axis  [16].  The  solution  for  the  region  x  >  a  along  x-  direction  is 
given  by 

X(x)  =  AeiK*x  +  Be~iK*x  (93) 

In  the  channel,  the  wave  moving  from  the  right  towards  negative  ^-direction  is  represented 
by  e~tKxX.  Since,  we  consider  the  case  of  the  wave  propagating  towards  positive  ^-direction 
only  (B  =  0),  we  have 

X(x)  =  AeiKxX  (94) 

The  above  equation  represents  a  traveling  wave  along  positive  for  x  »  a.  This  means  that 
the  propagation  wave  vector  has  to  be  positive,  that  is 

4  >  0  (95) 

For  the  region  x  »  a,  the  wave  function  in  the  transverse  direction  is  given  by 


Y(y)  =  Ccos(Kyy)  +  D  sin(nyy) 

Using  the  boundary  conditions  at  y  =  ±a  for  |x|  >>  a,  we  get 

Y(a)  =  0  =  C  cos (nya)  +  D  sin (Kya) 
Y(—a)  =  0  —  C  cos (nya)  —  D  sin (Kya) 

Adding  and  subtracting  above  equations,  we  get 

2  C  cos  (nya)  =  0 
2D  sin (nyd)  =  0 

Also 


co s(nya)  =  0 
nir 

ia  =  i 

nir 


n  =  1,  3,  5, 


(96) 


(97) 

(98) 


(99) 

(100) 


2  a 


(101) 
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Again  from  the  above  equation,  h/0  gives 


sin  (nya)  =  0 
nn 

Kva  =  n  =  0>2,4, 
nn 
2  a 


(102) 


For  n  =  0,  Y(x)  =  Const.,  which  is  a  trivial  solution.  Hence  for  both  odd  and  even  cases, 
we  have 


nn 


Ky=~,  n  =  1,2,  3, 4. 


Ky 


2  a 
nn 

w 


(103) 


where  w  =  2 a  is  the  width  of  the  channel.  The  free  propagation  along  the  channel  with 
transverse  quantization  is  a  typical  feature  of  quantum  wire.  The  value  of  the  transverse 
propagation  constant  is  minimum  when  n  —  1  and  is  given  by 

,  .  n 


\K. 


y  )min 


W 


(104) 


The  energy  corresponding  to  the  lowest  quantization  state  along  redirection  and  free  prop¬ 
agation  along  x-direction  is  given  by 

E  ~  EC  +  Kv)min ) 


2m 

K2 


E>—  Kt+  - 


.  n , 


o  1-*  .  V  2  ,  (105) 

2m  V  x  Kw  ) 

There  is  no  boundary  condition  that  restrict  the  values  of  kx.  ft  can  have  continuum  of  real 
values  greater  than  or  equal  to  zero.  Since  k2x  >  0,  the  lowest  energy  that  can  propagate  off 
to  infinity  along  the  channel  is  given  by  putting  k,x  =  0  in  the  above  expression. 

h2  n2 


Ef  = 


(106) 


2  mw2 

This  is  the  threshold  energy  for  the  propagation  along  the  channel  [1] .  The  threshold  energy 
does  not  depend  on  the  angle  between  the  crossed  wires.  This  is  obvious  because  inside  an 
infinite  channel,  which  is  far  from  the  crossed  intersection,  there  is  no  tunneling  from  the 
other  channel.  Any  state  with  energy  less  than  Et  will  give  the  bound  state.  Hence,  the 
condition  of  existence  of  localized  states  for  energy  E  within  the  channel  is 


E  <  Et 


(107) 
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In  our  approach,  the  energy  E  is  obtained  by  variational  method,  so  it  correspond  to  the 
minimum  energy  Ernm.  Many  literatures,  discussed  earlier,  show  that  for  a  single  straight 
wire  in  two  and  three  dimensions,  localization  is  possible  only  in  the  presence  of  corners 
or  the  bents  [9].  The  bents  and  corners  subsequently  lower  the  energy  of  the  system  due 
to  quantum  effects  and  hence  trapped  states  are  formed.  The  energy  corresponding  to 
the  bound  states  in  the  crossed  channels  is  the  bound  state  energy.  Also,  we  have  used 
variational  principle  to  estimate  the  energy,  hence  the  energy  is  the  ground  state  bound 
energy.  It  is  given  by 

Eb  =  Et  —  Emin  (108) 

The  choice  of  the  energy  scale  is  completely  arbitrary.  In  our  calculations,  for  convenience, 
the  origin  of  the  energy  scale  was  chosen  as  the  threshold  energy  Et.  In  such  a  case,  the 
bound  state  will  simply  be  given  by 

Eb  =  —Emin  (109) 

The  propagation  threshold  Et  is  a  constant  for  the  given  width  of  the  channels. 


B.  Calculation  of  the  bound  state  energy 


According  to  the  variational  principle,  the  expectation  value  of  the  Hamiltonian  is  mini¬ 


mized  with  respect  to  the  variational  a.  This  means  that  the  first  differential  of  the  expec¬ 


tation  value  of  the  Hamiltonian  (( H(a ))  )  with  respect  to  the  variational  parameter  (  a  ) 
must  vanish  at  some  value  a o  (say).  Hence,  from  Eqn.82,  we  get 


d(H(a)) 
da  |ao 

3 h2  d  /a2  +  2a  +  3\ 
mb2  da  \  11a  +  6  ) 

(llao  +  6)  (2ao  +  2)  —  (ajj  +  2ao  +  3)  11 
(lla0  +  6) 

llap  +  12a0  —  21 


0 

0 

0 

0 


a0  =  0.94001  (110) 


<92  (H(a)) 

da2 


l<*0 


>  0. 


Also,  we  find  that 
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This  cto  is  the  value  of  variational  parameter  which  minimizes  the  expectation  value  of 
Hamiltonian.  And  the  minimum  energy  is  obtained  by  substituting  ck0  in  {H),  which  is 
given  by 

E 'min 


The  threshold  energy  is  found  to  be 


Since  the  maximum  value  of  A  is  one,  Emin  is  less  than  the  threshold  energy  Et.  The  bound 
states  exist  in  the  classically  unbound  system  of  cross  wires  which  is  purely  due  to  quantum 
mechanical  effect.  This  minimum  energy  corresponds  to  the  binding  energy  of  an  electron 
in  the  cross  wires.  Hence,  the  bound  state  energy  (Eb)  of  the  electron  taking  the  threshold 
energy  as  zero  level  is  given  by 

Eb  =  —Ernin  =  — 1.058 -^-rA2 

rricE 

f)2 

=  —4.232 - -  cos2(di  —  6*2 ) 

mwz 

Ti2 

=  -4.232 - -cos  2  6  (113) 

mwz 

For  the  symmetrically  crossed  wire  system,  6  =  0  and  energy  Ernm  =  0.8575 Et.  Here,  Et  is 
the  unit  of  energy  representing  a  factor  equal  to  (n2 K2) / (2mw2)  (see  Eqn.112).  This  unit  is 
chosen  because  we  have  not  specified  the  nanowire  material  and  there  is  no  particular  values 
for  effective  mass  (m)  and  the  width  (w).  The  angular  dependence  of  the  bound  state  energy 
can  be  expressed  in  terms  of  the  threshold  energy  as 

Eb  =  -0.8575Etcos26  (114) 

The  typical  values  of  bound  state  energies  for  the  typical  rotation  angles  30°,  45°  and  60° 
are  —0.6431 Et,  —0A287Et  and  —  0.21437A  respectively. 
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C.  Angular  dependence  of  the  bound  state  energy 

In  our  calculations,  trial  wave  functions  were  so  chosen  that  they  vanish  outside  the  rhom¬ 
bus  region  even  in  the  channels.  The  variation  parameter  ’’a”  was  found  to  be  independent 
of  the  relative  orientation  of  the  channels.  This  means  that  ”a”  has  no  any  dependence  on 
the  angle  of  rotation  ( 0 ).  The  upper  bound  to  the  ground  state  energy  was  found  to  depend 
on  the  squared  cosine  of  the  angle  by  which  one  channel  is  rotated  with  respect  to  the  other. 
The  angle  of  rotation  0  represent  the  relative  rotation  (9\  —  62)  between  the  two  arms  of  the 
cross-wire  system  (see  figure  24).  When  0  =  0,  that  is,  if  each  arms  of  the  cross-wire  system 
is  rotated  by  equal  angle  with  respect  to  the  fixed  axis,  the  wires  intersect  at  right  angle. 
This  case  of  classically  unbound  system  is  known  to  posses  a  quantum  bound  state  [1]. 
This  very  fact  is  supplemented  by  the  model  (skewed  geometry)  formulated  in  this  thesis 
project.  It  is  clear  from  our  results  that  when  two  channels  intersect  exactly  at  right  angle 


FIG.  36.  Variation  of  bound  state  energy  of  an  electron  with  the  angle  of  rotation  8  =  8\  —  82- 
The  angle  of  rotation  ( 8 )  is  shown  in  the  adjoining  figure. 

(0  =  0),  the  system  is  in  strongly  bound  state  and  for  all  other  angles  the  system  is  less 
strongly  bound  (see  figure  36).  This  is  because,  in  the  absence  of  rotation  (0  =  0),  the  trial 
wave  function  for  an  electron  is  localized  well  inside  the  intersection  region  (see  figure  25). 
The  probability  of  finding  an  electron  is  maximum  inside  the  square  region  shown  by  sharp 
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kinks  in  the  density  plot  (see  figure  26).  In  the  case  of  finite  rotation  angle,  the  trial  wave 
functions  get  spread  out  of  the  intersection  region  due  to  change  in  the  shape  of  the  corners. 
This  feature  is  depicted  in  the  figures  [Fig.27-Fig.32],  When  two  channels  overlap  (6  =  7t/2), 
the  wave  function  is  confined  to  a  single  channel  (see  figure  33).  The  probability  densities 
becomes  uniform  in  the  channels.  The  probability  of  finding  an  electron  in  such  a  case  will 
be  anywhere  from  — oo  to  +oo  along  the  channel  (see  figure  34).  The  bound  state  energy 
is  zero  and  that  the  electron  can  never  be  trapped  in  a  single  straight  channel.  This  means 
that  an  electron  becomes  a  free  particle.  Thus,  bound  state  energy  of  an  electron  trapped 
in  a  cross  nanowires  varies  with  the  angular  deviation  ( 9\  — 92  =  9 )  as  an  integral  multiple 
of  7T.  The  bound  state  energy  corresponding  to  an  angle  9  of  rotation  is  identical  to  that  of 
rotation  by  an  angle  (n  —  9).  This  energy  is  symmetric  about  the  vertical  line  through  the 
angle  of  rotation  of  n/2.  This  must  be  true  because  the  geometry  of  cross- wire  system  for 
the  two  angles  differing  by  tt  is  exactly  the  same. 


FIG.  37.  Variation  of  bound  state  energy  with  the  angle  of  intersection  of  the  cross-wire  system 
u ;.  The  angle  oo  is  shown  in  the  adjoining  figure. 

The  dependence  of  the  bound  state  energy  with  the  angle  of  intersection  is  shown  in  the 
figure  37.  The  bound  state  energy  is  expressed  in  the  units  of  the  threshold  energy  ( Et )  for 
free  propagation.  It  is  clear  from  this  plot  that  the  bound  state  energy  varies  as  squared  sine 
of  the  intersection  angle  of  the  cross-wire  system.  This  energy  variation  is  symmetric  about 
oo  =  7t/2  and  odd  multiple  of  it.  The  variation  of  bound  state  energy  with  the  intersection 
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angle  is  complementary  to  that  of  energy  variation  with  the  angle  of  rotation  (see  figure 
36).  The  energy  corresponding  to  two  angles  of  intersection  differing  by  180°  is  identical. 
This  fact  is  obvious  from  the  geometrical  consideration.  The  geometry  for  the  intersection 
angle  u  is  exactly  the  same  as  its  supplementary  angle.  The  bound  state  energy  approaches 
maximum  value  as  the  angle  of  intersection  goes  to  90°.  Also,  as  the  origin  of  the  energy 
scale  is  taken  as  the  threshold  energy  itself,  bound  state  energy  is  zero  when  oj  —  0°.  This 
means  that  an  electron  can  freely  propagate  along  the  channels  of  the  cross-wire  system. 

Thus,  the  bound  state  energy  of  an  electron  trapped  at  the  skew-symmetric  cross-wire 
intersection  varies  as  squared  sine  of  the  angle  of  intersection  which  is  equivalent  to  the 
squared  cosine  of  the  relative  rotation  of  the  channels.  Analytically,  the  ground  state  bound 
energy  of  an  electron  for  our  model  system  can  be  written  as 

h 2  ft2 

Eb  =  -4.232 - -  cos2  6  =  -4.232 - -  sin2  u  (115) 

mw 1  mw2 

where,  the  rotation  angle  ( 6 )  and  the  intersection  angle  (cu)  are  related  by 

0  +  uj= |  (116) 

The  squared  sine  dependence  of  the  ground  state  bound  energy  on  the  cross-wire  intersection 
angle  implies  that  there  exist  unique  solution  for  the  angular  domain  [7t/2  >  co  >  0].  The 
system  considered  here  represent  both  bound  and  unbound  states.  The  existence  of  bound 
states  is  mainly  determined  by  the  angle  of  the  intersection  of  the  cross-wire  system.  This 
also  means  that  electron  is  trapped  by  the  presence  of  corners  and  bents.  If  corners  and 
bents  do  not  exist  as  in  the  case  of  a  single  wires  ( 9  =  vr/2),  then  the  system  is  unbounded. 
In  case  of  finite  rotation  angle  ( 6  ^  7t/2)  there  exist  bound  state  due  to  quantum  effects 
which  has  no  any  classical  counterpart. 
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VI.  RESULTS  AND  DISCUSSIONS 

A.  Summary  of  the  results 

The  present  work  can  be  summarized  as  follows: 

1.  We  have  found  that  the  bound  state  energy  of  an  electron  trapped  at  the  intersection  of 
the  quantum  nanowires  in  the  skewed  geometry  depend  on  the  angle  of  intersection  of 
the  channels.  The  ground  state  bound  energy  varies  as  squared  cosine  of  the  relative 
rotation  of  the  channels.  The  bound  state  energy  corresponding  to  an  angle  9  of 
rotation  is  identical  to  that  of  rotation  by  an  angle  (n  —  9).  In  fact,  it  must  be  true 
because  the  geometry  of  cross-wire  system  for  the  two  angles  differing  by  n  is  exactly 
the  same.  Also  parity  conservation  laws  require  that  the  energy  be  a  true  scalar,  hence 
it  must  be  invariant  under  reflection  symmetry. 

2.  The  results  of  our  calculations  show  that  the  role  of  tilted  nanowire  geometry  is  to  de¬ 
crease  the  bound  state  energy.  This  is  due  to  the  fact  that  in  the  tilted  configuration, 
one  pair  of  the  opposite  corners  become  sharper  than  90°  while  the  other  opposite  cor¬ 
ners  become  shallower,  making  the  system  elongated  in  one  direction,  thus  delocalizing 
the  electron.  In  the  absence  of  any  rotation  of  cross-wires  the  angle  of  intersection  is 
90°.  This  is  a  classically  unbound  system  and  is  known  to  have  quantum  bound  states 
due  to  the  presence  of  corners,  which  was  previously  reported  by  Schult  et  al.  [1].  Our 
calculations  show  that  for  any  finite  rotation  of  the  nanowire,  no  matter  however  small 
it  may  be,  there  is  always  a  decrease  of  the  bound  state  energy.  This  is  because  in  case 
of  no  rotation,  all  four  corners  are  identical  (included  angle  is  90°)  and  the  spatial  ex¬ 
tend  of  an  electron  wave  function  is  well  within  the  square  region  (see  figure  25).  But 
in  case  of  finite  rotation,  the  electron  wave  functions  get  spread  out  into  the  channels 
of  the  cross-wire  system  due  to  change  in  the  geometrical  constraints  (the  shape  of  the 
corners),  making  the  system  expanded  and  therefore  less  bound.  These  features  can  be 
visualized  directly  from  the  figures  [Fig. 27-Fig. 34]  in  Chapter  III.  Another  remarkable 
property  is  that,  when  the  angle  of  rotation  is  n/2  or  an  odd  integral  multiple  of  it, 
the  wave  function  is  further  spread  out  along  the  channel,  the  system  is  no  longer  a 
cross-wire  but  a  single  wire  (see  figures  33  and  34).  This  is  a  case  of  free  propagation 
and  the  quantum  bound  states  no  longer  exist  (see  figure  36).  The  transition  from 


quantum  bound  state  to  unbound  state  due  to  rotation  of  the  nanowire  is  the  intrinsic 
feature  which  can  be  useful  for  device  applications. 

3.  Different  geometries  of  the  cross-wire  systems  explored  earlier  (all  crossed  at  right 
angles)  show  the  effects  of  the  bound  states  on  the  electron  transport  properties  such 
as  tunneling,  transmission  and  reflection  [1,  7].  Various  ways  of  observing  the  bound 
states  were  also  suggested  through  measurement  of  the  conductance  using  two  dimen¬ 
sional  electron  gas  model  [2],  The  results  of  our  model  system  (cross- wire  intersecting 
at  an  arbitrary  angle)  not  only  supplement  the  theory  of  quantum  bound  state  in 
a  classically  unbound  system  but  give  the  quantitative  explanation  of  the  angular 
dependence  of  the  bound  state  energy. 

In  order  to  arrive  at  these  results,  many  steps  were  taken  in  sequential  pattern.  They  are 

•  Firstly,  after  extensive  study  of  the  literatures  on  the  cross-wire  systems  in  two  dimen¬ 
sion,  the  problem  that  takes  into  account  the  angular  dependence  of  the  bound  state 
energy  for  an  electron  was  determined  as  a  distinct  one. 

•  Secondly,  the  model  was  set  up  using  two  orthogonal  coordinate  systems  rigidly  fixed 
on  each  nanowires  and  orthogonal  transformation  was  applied  to  obtain  new  intersec¬ 
tion  geometry  (rhombus). 

•  Thirdly,  using  the  variational  principle  with  suitable  trial  wave  functions  for  the  oblique 
geometry,  calculations  were  carried  out  to  estimate  the  upper  bound  to  the  ground 
state  energy  for  an  electron  located  inside  the  cross-wire  system. 

•  Lastly,  the  results  were  visualized  using  Mathematica  programming  package. 


B.  Potential  applications 

As  it  has  been  stated  above,  our  results  supplement  a  theory  of  quantum  bound  states 
in  classically  unbound  systems  and  may  be  useful  to  interpret  electron  transport  peculiar¬ 
ities  in  realistic  systems  such  as  semiconductor  nanowire  films,  intersecting  nanorods  and 
carbon  nanotube  bundles  [7] .  Some  of  the  potential  applications  of  this  thesis  work  can  be 
summarized  as  follows: 
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1.  Transistors  and  switching  devices.  The  model  formulated  here  explains  the  effect  of 
intersection  angle  of  nanowires  on  the  bound  states  at  the  crossed  region.  The  bound 
state  energy  for  an  electron  at  the  crossed  intersection  decreases  as  the  angle  be¬ 
tween  the  nanowires  decreases.  This  feature  of  the  cross-wire  system  can  be  exploited 
to  design  ultra-sensitive  transistors  and  switching  devices  with  carbon  nanotubes,  in 
particular.  In  practice,  carbon  nanotube  bundles  can  be  regarded  as  a  collection  of 
asymmetrically  crossed  nanowires.  The  model  developed  here  could  be  used  to  inter¬ 
pret  electron  transport  peculiarities  in  such  nanotube  bundles.  Since  nanotubes  are 
very  rigid  structures  and  can  carry  high  densities  of  electric  current  without  losing 
their  properties,  the  crossed  nanotubes  structure  could  be  a  potential  candidate  for 
devices  such  as  single-electron  transistors  [19].  Recently,  resonances  have  been  ob¬ 
served  in  the  metallic  and  semiconducting  carbon  nanotubes  intersection  confirming 
the  existence  of  localized  states  [7] . 

2.  Amplifiers.  Semiconductor  quantum  wires  in  crossed  configuration  provide  superior 
performance  of  semiconductor  optical  amplifiers  (SOAs)  as  opposed  to  the  ampli¬ 
fiers  built  on  quantum  dots,  and  so  they  are  potential  candidates  for  future  solid-state 
quantum  information  and  communication  devices  [20].  With  the  development  of  litho¬ 
graphic  and  etching  techniques,  it  is  now  possible  to  manipulate  cross-wire  semicon¬ 
ductors  to  design  amplifiers.  The  model  suggested  here  could  be  of  immense  value  for 
controlling  the  quality  factor  of  such  amplifiers.  This  is  because  the  crossed  nanowires 
can  be  interchanged  from  bound  states  to  unbound  states  simply  by  varying  the  angle 
of  intersection.  This  property  of  the  cross-nanowire  system  could  be  used  to  develop 
power  efficient  amplifiers. 

3.  Photo-detectors.  The  theory  developed  here  can  be  used  to  study  the  performance  of 
quantum  wire  photodetectors.  In  the  crossed  nanowire  system,  an  electron  is  bound 
to  the  intersection  region  at  some  energy  level.  This  electron  can  be  excited  by  the 
light  radiation  (photon)  of  suitable  wavelength  that  scales  with  electron  bound  state 
energy.  The  photocurrent  will  be  produced  as  a  result  of  these  light  illumination  of  the 
system.  These  essential  properties  can  be  exploited  to  design  photo-detectors  which 
can  be  manipulated  and  controlled  by  varying  the  bound  state  energy  of  an  electron 
at  the  cross-wire  intersection. 
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4.  Lasers.  Lasers  may  be  designed  from  quantum  wires  which  have  desirable  properties 
with  regard  to  their  performances  [21].  Quantum  wire  lasers  are  predicted  to  have  a 
high  efficiency,  low  threshold  current  densities,  and  low  temperature  dependence  of 
the  threshold  current  [22].  ft  has  been  proved  that  quantum  wire  lasers  have  lower 
threshold  current  density  as  compared  to  the  quantum  dots  lasers  [23].  The  energy 
spectrum  of  the  crossed  quantum  wires  can  be  varied  over  a  wide  range  by  simply 
changing  the  angle  of  intersection  between  them.  The  population  inversion  can  be 
achieved  by  controlling  the  voltage  across  the  crossed  wires  and  also  varying  the  in¬ 
tersection  angle.  Thus,  the  use  of  crossed  quantum  wires  may  offer  possibilities  for 
power-efficient  nonlinear  devices  in  semiconductor  technology. 

5.  Solar  cells.  The  theory  developed  here  can  be  used  to  design  the  energy  conversion 
devices  such  as  solar  cells.  Quantum  wire  solar  cells  have  the  potential  to  increase 
the  maximum  attainable  conversion  efficiency  of  solar  photon  conversion  by  utilizing 
hot  photogenerated  carriers  to  produce  higher  photovoltages  or  higher  photocurrents. 
The  nanowire  geometry  provides  potential  advantages  over  planar  wafer-based  or  thin- 
film  solar  cells  in  every  step  of  the  photoconversion  process  [24],  Knowledge  of  the 
quantum  bound  states  in  the  cross-nanowire  system  is  useful  to  explain  the  conversion 
of  photon  flux  into  an  electron  flux  through  the  nanowire  channels.  In  the  crossed 
geometry,  the  nanowires  will  have  greater  resistance  to  the  flow  of  electron  than  the 
corresponding  parallel  configuration.  Therefore,  understanding  the  properties  of  the 
crossed  nanowire  system  would  allow  one  to  improve  the  efficiency  of  the  quantum 
wire  solar  cells  as  compared  to  the  conventional  solar  cells. 

Finally,  it  can  be  concluded  that  the  existence  of  bound  states  in  the  cross-wire  system, 
which  is  classically  unbound  system  demonstrates  that  the  laws  governing  microscopic  world 
are  inherently  quantum  mechanical.  Extension  of  this  work  to  incorporate  more  realistic 
models  of  the  confining  potential,  inhomogeneous  external  fields  and  temperature  fluctua¬ 
tions  appear  fascinating  but  still  remains  to  be  done. 
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VIII.  APPENDIX  I 


A.  Calculation  of  the  normalization  constant 
For  oblique  quadrants  Rn  For  the  oblique  quadrants  defined  by  Rn  where,  n  =  1,4, 


the  trial  functions  are: 


|  x'y'1 


$(x',y')  =  A  (  1 - —  ]  e  °  for  \x'\  <  b,  \y'\  <  b 


(117) 


For  an  oblique  quadrant  Ri.  We  take  the  trial  function  as 


$1(x',y>)=A(l-'^f)  e~a 


(H8) 


and  the  region  of  integration  are  from  x'  =  0  to  x'  =  b  and  if  —  0  to  y'  —  b.  Hence, 


Ri=  /  &l(x',y')  &i(x',y')  secOdx'dy 

J x'=0  J y'= 0 

pb  pb 

=  A2  seed  / 

J x'=0  J y'= 0 

rb  rb 


2  J  "  1  1  ‘l-XJ)\-2*dx'dy' 


/  /  \  2 

=  A2  sec  6  e~2a  I  I  (  1  —  I  dx’  dy’ 


J  x'=0  J  y'= 0 
pb  pb 

=  A2  sec  9  e~2a  /  /  I  1 

J x'=0  J y'= 0 


b2 


2  x'y'  x,2yl2\  ,  , 


b2 


+ 


b 4 


dx'  dy' 


=  A2  sec  9  e~2a  b2  - 


2  b2b2  b3  b3 
4  b2  +  9  V 


7,2  1 2 

=  A2  sec  6  e~2a  (  b2  -  —  +  — 

2  9 


=  A2b2  sec  6  e~2a  (-  +  - 
,2  9 


=  A2  If2  sec  9  e~2a  (  — 


18 


11 


=  —  A2  b2  sec  6  e~2a 
18 


For  an  oblique  quadrant  R-2:  We  have  trial  function  as 


4>2(z',j/)=7l  (1  +  ^)  e-“ 


(119) 


(120) 


and  the  region  of  integration  are  from  x'  =  —  b  to  x'  =  0  and  y'  —  0  to  y'  —  b.  Hence, 

r0  rb 

R-2  —  I  /  ®*2(x',y')  $2(x',y')  seeddx’  dy’ 

J x'——b  J  y'= 0 
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=  A2  seed  [  [  +  e~2a  dx'  dy' 

J x'=—b  Jy’= 0  V  V  J 

=  —A2  secde~2a  [  [  (l  +  Kf)  dx' dy' 

Jx'= 0  Jy’= 0  V  b  J 

Let  u  =  —x',  so  that  du  =  —dx'  ,  therefore,  when  x'  =  0,w  =  0  and  when  x'  =  —  b,  u  =  b. 
Therefore, 

R-2  =  A2  sec  6  e~2a  [  [  (l  —  dudy'  —  R\  (121) 

Ju= o  Jy'=  o  V  by 

For  an  oblique  quadrant  R,:> :  We  take  the  trial  function  as 

^3  (x',y)  =  A  ^1  -  exp  (—a)  (122) 

and  the  region  of  integration  are  from  x'  =  —  b  to  x'  =  0  and  y'  =  —b  to  y'  =  0.  Hence, 

R'i  =  [  [  mx',y')$3(x',y')  sec  Odx'dy' 

J  x'=—b  J  y'=—b 

=  A2  seed  [  [  ( 1  -  e~2adx'dy' 

J  x'=—b  J  y’=—b  \  b  J 

=  A2  sec  6  e~2a  f  [  (l  -  dx' dy' 

Jx'= 0  Jy'= o  V  by 

Let  u  =  —x',  so  that  du  =  —  dx'  ,  therefore,  when  x'  =  0,  u  =  0  and  when  x'  =  —6,  u  =  b. 
also  let  v  =  —y\  so  that  dv  =  —  dy'  ,  therefore,  when  y'  =  0,  v  =  0  and  when  y'  =  —b,  v  =  b. 
Therefore, 

R3  =  A2  sec  6  e~2a  [  [  (l — — d/a  dv  —  R i  (123) 

Ju= o  Jv= o  1  b  J 

and  hnally,  For  an  oblique  quadrant  /i'4 :  We  have  the  trial  function  as 

<h4 {x',y')  =  A  ^1  +  e_"  for  x'  >  0,  y'  <  0  (124) 

and  the  region  of  integration  are  from  x'  =  0  to  x'  =  b  and  y'  =  —b  to  y'  =  0.  Hence, 

nb 

R4  —  I  /  $\{x',  y')  $4(x',  y')  sec  6  dx'  dy' 

J x'=0  J  y'=—b 

=  A2  seed  [  [  (l  +  ^Tf)  e~2adx'dy' 

Jx'= 0  J  y'=—b  \  by 

=  -A2  sec  d  e~2a  f  [  +  dx' dy' 

Jx'= 0  Jy'= 0  V  by 


65 


Let  v  =  —y',  so  that  dv  =  —dy'  ,  therefore,  when  y'  —  0,u  =  0  and  when  y'  =  —  b,  v  =  b. 
Therefore, 

i?4  =  A2  sec  9  e~2a  f  f  (^1  —  — ^  dx'  dv  —  R\  (125) 

Jx'= o  Jv= o  V  b2  ) 

The  total  contribution  to  the  normalization  integral  from  the  regions  Rn  is  given  by 

4 

i=i 

=  4i?i 

=  4  A2b 2  sec  9  e~2a 

=  2^-A2b2  secde~2a  (126) 

9 

For  the  open  channels  Cn  For  the  open  channels  defined  by  Cn  where,  n  =  (1,  8) ,  the 
trial  functions  are: 

4>(x',y')  =  A  ^1  —  for  |a/|  >  b ,  \y'\  <  b 

=  A  ^1  —  e~~^  for  \y'\  >  b ,  \x'\  <  b  (127) 

For  the  channel  C\ :  We  get 


Mx',y')  =  A^  1-|)  e~^  (128) 

and  the  region  of  integration  are  from  xf  =  b  to  x'  =  oo  and  y'  —  0  to  y'  —  b.  Hence, 

/»oo  rb 

Ci—  /  fi(x',  y')  <$>i(x' ,  y')  sec  6  dx1  dy' 

J x'=b  J y'= 0 

=  A2  seed  [  [  fl-f)  e-^dx'dy' 

Jx'=b  J y'=0  \  b  J 


'  r°°  „  / 

rh  f  2 

A2  seed 

[  e~  V  dx' 

/  f1-  7)  dy' 

_Jx'=b 

Jy'=0  V  b  J 

—  (0  -  e~2a) 
2  a  v  ’ 


=  A2  sec  6 

=  A2  seed  —e~2a 
2  a 

A2 

=  —  b2  sec  9  e~2a 
6  a 


For  the  channel  Cx.  We  have 


(129) 


<h(x',y') 


(130) 
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and  the  region  of  integration  are  from  x'  =  0  to  x'  =  b  and  y'  —  b  to  y'  —  oo.  Hence, 

pb  poo 

C-2  —  I  /  (j)l(x',y')  4>2{x\y')  secddx'dy' 

J x'=0  J y'=b 


=  A2  sec  9 


=  A2  sec  9 


I  x'=0  J  y'=b 


'  /  \  2 

x  \  _W  ,  ,  , 
1  — —  J  e  b  ax  ay 


f--j  dx’ 


2  a  y  , 

e  dy 


=  A2  sec  9  -  —e  2a  —  C\ 
3  2a 


For  the  channel  Ch:  We  have 


03  (x,y)  =  A\[1  +  Jje  b 

and  the  region  of  integration  are  from  x’  =  —  b  to  x1  =  0  and  y'  =  b  to  y'  =  oo. 


Hence, 


p  0  poo 

C3=  /  cj)*3(x\y')(j)3(x\y')  seed  dx' dy' 

J  x'=—b  J  y'—b 


—  A2  sec  9 


r»0  poo 


'  x'=—b  J  y'=b 


'  /  \  2 

X  \  2a  y' 

1  +  —  I  e  S;  dx  d/y 


=  —A2  sec  6 


i  +  -J  dx1 


2 a  yr  , 

e  b~  dy 


(131) 


(132) 


(133) 


Let  u  =  —x',  so  that  du  =  —dx'  ,  therefore,  when  x'  =  0,u  =  0  and  when  x'  =  —  b,  u  =  6. 
Therefore, 


C3  =  H2  sec  6 


1  -  !)2  du 


2a  yr  , 

e  ~b  -  dy 


=  A2  seed  [-1  [—  e-2“l  =  Ci 

_3j  [2a 

For  the  channel  C4:The  trial  function  is  taken  as 


(134) 


M* .  v )  =  A  ( i  -  j  j e  * 


(135) 


and  the  region  of  integration  are  from  x'  =  — oo  to  x'  =  —b  and  y’  —  0  to  if  —  b.  Hence, 


p—b  pb 

C4  =  /  /  04  (x',  y')  04(a)',  y')  sec  d  dx'  dy' 

J  x'=—oo  Jy'= 0 
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=  A2sec9  [  [  (\  e2¥'rdx'  dy' 

Jx'=- oo  Jy'= 0  V  °  J 


A2  sec  9 

[  dx' 

f  dy' 

_Jx'——b 

Jy'= 0  V  b  J 

Let  u  =  —x',  so  that  du  =  —dx'  ,  therefore,  when  x'  =  —b,  u  =  a  and  when  x'  =  —  oo 
Therefore, 


Ca  =  A2  sec  9 


=  A2  sec  0 


—  2a.  u 

e  b  du 


'y'= o 


-2a 


2a 


=  C, 


For  the  channel  CV :  The  trial  function  is  taken  as 


2/0  =  A  I  1  +  \  )  e“fc 


and  the  region  of  integration  are  from  x'  =  —  oo  to  x'  =  —b  and  y'  =  —b  to  y'  =  0. 

p—b  /*0 


cv  = 


05  (af ,  y  )  05  (x1 ,  y  )  sec  9  dx' dy' 


'  x'=— oo  J  y'=—b 

r—b  rO 


=  A2  sec  0 


=  A2  sec  9 


x'=—oo  Jy'=—b 
r-b 


1  +  —  e  &  dx  dy 
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2ol  x1  . 

e  b  dx 


'  X—  —  OC 


1  +  y  )  dy' 


Let  u  =  —x',  so  that  du  =  —dx'  ,  therefore,  when  a;'  =  —  oo,  u  —  oo  and  when  x'  = 
also,  let  v  =  —y',  so  that  dv  =  —  dy',  therefore,  when  y'  —  0,  v  —  0  and  when  y'  = 
Therefore, 


CV,  =  A2  sec  $ 


=  A2  sec  9 


—  2  a.  u 

e  b  du 


L  J  u=b 


-2  a 


2a 


L  J  t;=0 

=  C1 


For  the  channel  CV:  The  trial  function  is  taken  as 


(136) 
,  u  —  oo. 

(137) 

(138) 
ffence, 

(139) 

b,  u  =  b. 
■b,  v  =  b. 

(140) 


'klx'.  y')  =  A  (l  +  j 


(141) 


and  the  region  of  integration  are  from  x'  =  —  b  to  x'  =  0  and  y'  =  —  oo  to  y'  =  —b.  Hence, 
C6=  [  I  4>l(x',y')  (j)&(x',y')  secOdx'dy' 

J  x'=—b  J  y'=— oo 

=  A2  seed  [°  I*  b  ^ 

J  x'=—b  J  y'=— oo 


=  A2  sec  6 


r-'  =  -h  .lii'  —  - 

fO 


1  +  —  e  b  dx  dy 
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1  H - )  dx' 
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c—b 


2 cl  y'  . 

e  b  dy 


(142) 


Let  u  =  —x',  so  that  du  =  —dx'  ,  therefore,  when  x'  =  —b,  u  =  b  and  when  x'  —  0,  u  —  0. 
also,  let  v  =  —y'i  so  that  dv  =  —dy'  ,  therefore,  when  y'  =  —oo,v  =  oo  and  when  y'  = 
—b,  v  =  b.  Therefore, 


=  A2  sec  6 


=  A2  sec  9 


L  J  u= o 

b 
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f  ^Ydu 


_ 2a  v 

e  b  dv 


'  v=b 


-2a 


2  a 


=  Ci 


(143) 


For  the  channel  CV:  The  trial  function  is  taken  as 


x 


=  A  (  1  -  -  )  e 


(144) 


and  the  region  of  integration  are  from  x'  =  0  to  x'  =  b  and  y’  =  —  oo  to  y'  =  —b.  Hence, 

-fe  r-b 
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x'=0  J  y'=— oo 
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e  b  dy 
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(145) 


Let  v  =  —y',  so  that  dv  =  —dy’  ,  therefore,  when  y'  =  —b,  v  —  b  and  when  y'  =  —oo,  v  =  oo. 
Therefore, 


C7  =  A2  sec  9 


=  A2  sec  9 


’  x'=0 
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1  — —  1  dx' 
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—  2a  v 

e  b  dv 
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2  a 


—  C\ 


(146) 


For  the  channel  Cs:  The  trial  function  is  taken  as 


y 


x',y')  =  a  (  1  +  —  )  e 


(147) 


69 


and  the  region  of  integration  are  from  x'  =  b  to  x'  =  oo  and  y'  =  — b  to  if  =  0.  Hence, 

/»oo  /*0 

/  (j)*8{x',  y')  <p8(x',  y')  sec  9  dx’  dy 


C8  = 


1  x'=b  Jy'=-b 


=  A2  sec  9 


=  A2  sec  9 
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f  x'=b  J  y'=—b 


1  +  e  b  dx'  dy' 


—  2a.  x  . 

e  b  dx 


'  x'=b 


-0 


'y'=-b 


1  +  j  I  <l,J 


(148) 


Let  v  =  —y',  so  that  dv  =  —dy'  ,  therefore,  when  y'  =  0,v  =  0  and  when  y'  =  —b,v  =  b. 
Therefore, 


C8  =  A2  sec  9 
=  A2  sec  9 


—  2a  x  j 

e  b  dx 


L  J  x'=b 

b 


-2a 


2a 


U 

~  C\ 


1  —  -  )  dv 
b 


(149) 


The  total  contribution  to  the  normalization  integral  from  the  channels  Cn  is  given  by 


i= 1 

=  8  Ci 

'  A2 

=  8  [  —  b2  sec  9  e~2a 
6a 


4  A2 
3  a 


b2  sec  9  e  2a 


Hence,  from  the  normalization  condition  we  have, 


(150) 


I\  +  J\  —  1 

99  A  A2 

—  A2  b2  sec  9  e~2a  +  —  b2  sec  9  e~2a  =  1 
9  3a 


2  A2  b2  sec  9  e~2a 
9  a 


(11a +  6)  =  1 


H2  = 


9  a  e 


2  a 


2b2  (11a  +  6)  seed 


A  = 


9  a  e 


2a 


2  b2  (11a  +  6)  sec  9 


(151) 


This  is  the  required  normalization  constant  for  the  given  set  of  trial  functions. 
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IX.  APPENDIX  II 


A.  Calculation  of  channel  integrals 
For  the  regions  Cn  (n=l  to  8) 

For  channel  C\\  x'  =  b  to  x'  =  oo  and  y'  —  0  to  y'  —  b ,  x',  y'  >  0 


and 


TT  - h 2  /  d2  d2  n  d2 

H  =  - —  |  +  ttxx  -  2  sin  6 


2 m  V  dx'2  d y’2 


dx’dy 1 


(152) 


gms'K'lf 

dV3 

dx 12 

d2<fs  Aa  —ay' 
dx’dy’  b2 


=  0 


3  = 


—h2Aa 
2 m  b2 


x 


a  (  1  H - 1+2  sin  d 

b 


—  ay 

e  b 


Hence,  the  matrix  element  of  H  is  given  by 


Jci  =  ((pi(x'  ^y^lHlip^x'  ,y')) 

-h2A2a  fb  /  y'  .  _ 

2m62  Jx>=b  Jy,=0  V  & 


a  (  1  —  )  —  2  sin  d 

0 


'  x'=b 


-K2A2a 
2  mb2 

— K2A2a  sec  6 
2  mb2 

— K2A2e~2a  seed  r 
4m 


~a3:  ,  / 

e  &  ax 


r> 

ly'= 0  \  b 


a  (  1  — -  ]  —  2  sin  0 
b 


sec  Odx'dy' 
dy' 


"  b  e-2a 

_2a 

a|^-2smoH 


a  •  a 

- sm  t) 

3 


(153) 


For  channel  C2:  x'  =  0  to  x'  =  b  and  y'  =  b  to  y'  =  oo,  x' ,  y'  >  0 


,  ,  x'\  zW 

7^2  =  H.  (  1  —  J  6  b 


(154) 
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From  symmetry,  it  is  clear  that  Jc2  =  Jci- 

For  channel  C3:  x'  =  —  b  to  x'  =  0  and  y'  —  b  to  y'  —  oo,  x'  <  0,  y'  >  b 


.  .  „  X  \  ^v_ 

•fa  =  A  (  1  +  -  )  e 


(155) 


and 


s-yyiK1' ■“ 


d2<p3  A  (a\2  A  ,  y' 


=  A  -  1  +  —  e~b 

dy'2  \b)  [  b  1 

d2<p3 


dx'2 

<92F3 


=  0 


j\.0i  —oty' 

dx'dy'  b2 


3  = 


— h2Aa 
2  mb2 


x 


a  (  1  +  —  )  +  2  sin  6 
b 


—  ay 

e  b 


Hence,  the  matrix  element  of  H  is  given  by 

Jcs  =  {<f3(x',y')\H\<f3(x',y')) 


X 


a  (  1  +  —  ]  +  2  sin  6 
b 


x 


a  (  1  +  —  ]  +  2  sin  9 
b 


sec  6  dx'dy' 
dx' 


Let  u  =  —x',  so  that  du  =  —  dx'  ,  therefore,  when  x'  =  0,  u  =  0  and  when  x'  =  —6,  u  =  b. 
Therefore 


Jc3  — 


h2  A2  a  sec  6 
2  mb2 

—  h2A2asec6 

2  mb2 

—  h2A2asec6 

2  mb2 


-oo  „  ,  rO 

~2ay  ,  /  / 
e  f ’dy 


y'=b 
oo 


a  ( 1  —  —  )  +  2  sin  6 
b 


du 


—  2  <y.yr  , 

e  b  dy 


'  u— 0 


u 


a  ( 1  —  —  )  +  2  sin  9 
b 


du 


-2  a 


2  a 


a^)+2sin0(^ 


—k2A2e  2asec  9  \ 
Am 


a  ■  n 

— h  sm  t) 
3 


(156) 


For  channel  C4:  x'  =  —  00  to  x'  =  —6  and  y'  —  0  to  y'  —  b,  x'  <  0,  y'  >  0 


</?4 


(157) 


From  symmetry,  it  is  clear  that  Jca  =  Jcs- 

For  channel  C5:  x'  =  —  oo  to  x'  =  — b  and  y'  =  —b  to  y'  =  0,  x',  y'  <  0 
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A 


/ 

ax 

e  b 


and 


d<p5 

dx' 

dx'2 

d2y5 

d  y'2 
d2(f5 
dx'dy' 

H(p  5 


Aot  ai1 
- P  b 

b2 


—h2Aa 


2  mb2 


a 


2  sin  0 


Hence,  the  matrix  element  of  H  is  given  by 


(158) 


Jcb  =  {^{x\y')\H\y}b(x\y')) 

—h2A2a  rb  r° 


2  mb2 
-K2A2a  sec  6  f~b 


V  \  olx' 
1+  6  >e‘ 


a  (  1  +  7  )  —  2  sin  0 
b 


sec  Odx'dy' 


2  mb2 


2  ax'  . 

e  b  dx 


0  /1  y' 

1  +  T 


'  y  =— 


y 


a  (  1  +  —  I  —  2  sin  9 
b 


dy1 


Let  u  =  — x so  that  du  =  —  dx'  ,  therefore,  when  x'  =  —6,  u  —  b  and  when  x'  =  —  oo,  u  =  oo. 
also,  let  v  =  —y',  so  that  dv  =  —dy’ ,  therefore,  when  y'  =  —b,  v  —  b  and  when  y'  —  0,  v  —  0. 
Therefore, 


Jc5  — 


-h2A2a  sec  9  fb 


2  mb2 

— h2A2a  sec  9 
2  mb2 

— h2A2a  sec  9 


e  b  du 


'  U= OO 
poo 


—  2  OLU 

e  b  du 


2  mb2 


u=b 

■  b 
2  a 


-2a 


v=0 

b 

OL  I  - 


x-\ 


a  ( 1  —  -  )  —  2  sin  9 
b 


a  ( 1  —  -  )  —  2  sin  9 
b 


dv 

dv 


2  sin  9  - 
2 


— K2A2e  2asec  9  \ 


4  m 


a  ■  a 
- sm  9 

3 


(159) 


For  channel  C6:  x'  =  —  b  to  x'  =  0  and  y'  =  —oo  to  y'  =  — b ,  x'  <  0,  y'  <  0 


.  .  x'\  ay' 

7,6  =  24(1  +  —  J  e  b 


(160) 
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From  symmetry,  it  is  clear  that  Jc%  =  Jc5- 

For  channel  Cf.  x'  =  0  to  x'  =  b  and  y'  =  —  oo  to  y'  =  —b,  x'  >  0,  y'  <  0 

.A  x'\ 

w  =  A  V1  “  ~b )  e  b 


and 


a 


=  a  iFi2Ci-7W 


&P 7 
dy' 

92(P7 
dy'2 

dy7 

dx'2 
d2(P7 

dx'dy'  b 2 
-  —h2Aa 

^  ~  2 mb2 


x 


=  A  (  -  )  (  1  -  —  left 
b 

a 

b 


=  0 

Aoi  ay' 

=  — rre  b 


X 


a  (  1 - ]  +  2  sin  9 

b 


ay 

e  b 


Hence  the  matrix  element  of  H  is  given  by 


(161) 


Jc7  =  (f7(x',y')\H\ip7(x',y')) 

—h2A2a  rb  rb 
2 mb2  Jx,=0  Jy ,=.00 

— h2A2a  sec  6  f~b 


x'\  2a  y' 

1- —) e b 

b 


a  (  1  — -  )  +  2  sin  9 
b 


sec  ddx'dy' 


2  mb2 


,  / 

e  b  dy 


I  x'=0 


X 


a  |  1  — —  )  +  2  sin  9 
b 


dx' 


Let  v  =  —y',  so  that  dv  =  —dy'  ,  therefore,  when  y'  =  — oo,  v  —  oo  and  when  y'  =  —6,  v  =  b. 
Therefore, 


Jc7  — 


h2A2a  sec  9  fh 


2  mb2 

—  h2A2asec6 

2  mb2 

—  h2A2asec9 

2  mb2 


e  b  dv 


v=oo 

OO 


b  / 


—  2o:n 

e  &  cfo 


v=b 

'  b_ 
‘la 


'  tc'=0 

t 

'  x'=0 


X 


a  (  1  — —  )  +  2  sin  9 
b 


l~-b 


x 


a  |  1  — —  )  +  2  sin  9 
b 


dx ' 
dx' 


-2a 


a  x  +  2  sin  9  i  - 


-h2A2e  2asec9 
Am 


a  ■  n 

13+Sm$. 


For  channel  Cg-  x’  =  b  to  x1  =  oo  and  y'  =  —6  to  y'  =  0,  x1  >  0,  y'  <  0 

^(1  +  i)e"¥ 

From  symmetry,  it  is  clear  that  Jcs  —  Jc7  and  therefore 

Jci  —  Jc2  =  JC5  ~  Jc6 


(162) 


(163) 
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and 


Jc3  —  Jca  —  Jc7  —  Jcs 


The  total  the  matrix  element  of  H  is  given  by  for  the  channels  Cn  is  given  by 

(H)c  =  Jc\  +  Jc2  +  Jc3  +  JcA  +  Jc5  +  Jce  +  J C7  +  Jc8 

=  2  JCi  +  27^3  +  2Jc5  +  2Jc<7 


=  4  (Jci  +  Jc3) 

^  H2A2e~2a  sec  0  r 


4m 

2h2A2ae~2a  sec  9 
3  m 


a  ■  a 
- sm  9 

3 


+ 


— K2A2e  2a  sec  9  r 


4  m 


a  ■  n 

— h  sm  9 
3 


(164) 


B.  Calculation  of  rhombus  integrals 


For  the  regions  Rn  (n=l  to  4) 


For  oblique  quadrant  R\\  x'  —  0  to  x'  —  b  and  if  —  0  to  y'  —  b  x',  y'  >  0 


and  the  Hamiltonian  operator  is  given  by 

-h2  (  d2  d 2 


So  we  have 


H  = 


+ 


2m  \  dx'2  d y'2 

d(pi 


d2  \ 
2  Smddx'dy' ) 


=  Av'e-a 


dx'  b2 

aVi  A 

=  --e 


dx'  dy'  b2 
d2p\ 
dx'2 
d2p 


=  0 


dy 


J2 


l  =  0 


(165) 


(166) 


A  K2  sin  9 


m 


b2 
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Hence  the  matrix  element  of  H  is  given  by 

Iri  =  {<PiW,y')\H\<Pi(x',y')) 

Ah2  sind  „  fb  i 


m  b 2  J 

A2  h 2  sin  9  sec  9 
m  b 2 

A2  h 2  sin  9  sec  9 


A  |  1  _  ®  y  \  a - n  J 


m  b2 
3  h2  A2 


x'=0  J  y'= 0 

rb 

-2a  / 


’  x'=0  Jy'= 0 
-2a  I  ft 2 


b2 


1  - 


e  “  sec  9  dx  dy 


x'y' 

b2 


dx  dy' 


4  m 


sin  9  sec  9  e 


-2a 


For  oblique  quadrant  R-2'.  x'  =  —6  to  x'  =  0  and  y'  =  0  to  y'  =  b  x'  <  0,  y'  >  0 


Hence,  we  get 


<P2  =  A[l  +  ?f]  e- 


<9^2  71  ,  _c 


=  e 


dx'  b2 
A 


d2y> 


2  =  ^  e“Q 


Hlp2  = 


dx'  dy'  b2 
d2<f2 
dx'2 
d2(p2 
dy 12 

-h2 


=  0 
=  0 


2  m 
A  h2  sin  9 
m  b2 


2  sind(4e"a) 
b 2 


Hence,  the  matrix  element  of  H  is  given  by 


(167) 


(168) 


Ir2  =  {^2(x',y')\H\^2(x',y')) 


A  h 2  sin  9 
m  b2 
A2  h2  sin  9  sec  9 
m  b2 


x'=—b  Jy'=0 
fO 

-2  a  / 


A  1  + 


x'y' 


1  + 


e  a  sec  9  dx'  dy' 


x'y' 

~W 


dx'  dy' 


J  x'=—b  Jy'= 0 

Let  u  =  —x',  so  that  du  =  —dx',  therefore,  when  x!  =  0,u  =  0  and  when  x'  =  —  b,  u  =  b. 
Therefore 


Ir2  — 


A2  h2  sin  9  sec  9 
m  b2 


-2a 


u=b  Jy'=0 


b  /  mir 

1  -  —  )  dudy’ 
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=  —If 


(169) 


For  oblique  quadrant  i?3:  x1  =  —6  to  xf  =  0  and  y1  —  —  b  to  y1  —  0  x'  <  0,  y'  <  0 

^  =  x(l-^)«-  (170) 


<9^3  71  ,  _Q 

^  =  ~Wye 
dVs  =  _A  _a 

dx'  dy'  b 2  6 


A  H2  sin  6 


Hence  the  matrix  element  of  H  is  given  by 


Im  =  (ipz{x'  ,y')\H\y>z(x'  ,y')) 
Ah2  sin 9 


0  r  0 


'  x'=—b  J  y'=—b 


.All - —  J  e  a  sec  $  dx  dy 


=  -  ^  3il‘f  secg  e-2"  /'”  /”  (l  -  77)  dx’  dy’ 

m  b2  Jx,=_b  Jy,=_b  V  b2  J 

Let  u  =  —a/,  so  that  du  =  —  dx'  ,  therefore,  when  a/  =  0,n  =  0  and  when  x'  =  — b,u  =  b. 
also,  let  v  =  —y',  so  that  dv  =  —  dy1,  therefore,  when  y'  —  0,  v  —  0  and  when  y'  =  —6,  v  =  b. 
Therefore,  we  get 

A2h2  sind  sec  9  ,  , 


Iro  =  —  - 


sin  9  sec  0  2  /'°  /'°  /  uv  \ 

^ — ■  LL 


H2h2  sind seed  _2a  fb  fb 

™  b2  Ju= o  Jv= o 


1  -  du  dv 
b2  J 


(171) 


For  oblique  quadrant  iTy  a/  =  0  to  a;7  =  b  and  y'  =  —  6  to  7/  =  0  a;7  >  0,  y7  <  0 

,4^(l  +  f) 


(172) 
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dV±  =  -v'e-a 

dx'  b 2  U 
92  f  4:  =  A  ^ 


dx'  dy'  b 2 
92<Pa 
dx'2 
d2if2 


=  0 


dy 


/2 


=  0 


^4=  i(-2  sin9(^e‘“: 

A  h 2  sin  0 


m 


b2 


Hence,  the  matrix  element  of  H  is  given  by 


Ira  =  {va(x'  ,y')\H\yA(x'  ,y')) 


A  h2  sin  9 


-^e  *.  L  l  . A  0 + l~w  •  e_“  secedx'dy' 


A2  h2  sin#  sec  9  2  fb  f°  (  x'y'\  ,  , 

m*  e  LL>  O  + 


Let  v  =  —y',  so  that  dv  =  —dy',  therefore,  when  y'  —  0,  v  —  0  and  when  y' 
Therefore,  we  have 


Ira  — 


A2  h 2  sin  9  sec  6 
m  b2 

A2  h2  sin  9  sec  9 


-2a 


lx'= 0  J  v=b 


x'v 


—2a 


m  b2 


'x'=0  J  v= 0 


1  — —  I  dx1  dv 
b 2 


l-$]  dx'dv 
b 1 


=  -I 


Ri 


The  total  the  matrix  element  of  H  is  given  by  for  the  region  Rn  is  given  by 


(H)r  ~  Iri  +  Ir2  +  Irs  +  I Ra 
=  Iri  ~  Ira  +  Ir,  —  Ir, 


b,  v  =  b. 


(173) 


0 


(174) 
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X.  APPENDIX  III 


A.  Calculation  of  matrix  elements  for  outer  four  arms 


The  outer  four  arms  are  defined  by  the  following  conditions: 


x1  =  b 

to 

X1  = 

oo  ,y' 

=  0 

x'  =  —oo 

to 

x'  = 

- b,y ' 

=  0 

y’  =  b 

to 

y'  = 

OO,  X1 

=  0 

V  =  -oo 

to 

y'  = 

—b,  x1 

=  0 

For  the  right  arm,  x'  =  b  to  x'  =  oo,  y'  =  0  and  the  trial  function  near  its  proximity  is  given 
by 

/  b/l\ 

(175) 


I  y' 1 

4  =  A  (  1  - -j-  )  e  I 


The  Hamiltonian  operator  is  given  by 


„  -h2  (  d2  d 2  0  .  .  d 2 

H  =  - —  I  ^7  +  -  2  sin  9- 


(176) 


2m  \  dx'2  dy 12  dx'dy'  / 

As  there  is  no  discontinuity  in  the  derivative  of  ip\  with  respect  to  x' ,  contribution  is 
not  taken  in  the  calculation  of  matrix  element  of  H  in  this  roof  line  as  it  is  already  taken 
into  account  in  the  bulk  integration.  Hence,  we  have  to  consider  H  as 


2m  \dy'2  ’  dx'dy' 

Using, \y' \  =y'  [%')  -  9(-y%  and  5(y')  =  S(-y')  we  get, 


(177) 


4  44)  =  w  [%')+*(-»')] 


dy 


A  (\  -  Ml  e=¥  = 

dx'  \  b  /  b 


b 


Hip  i  =  A 


-h2 

2m 


—2  n 

T%')-  2  sin 9^  {9{y')  —  9(—y')) 


e  b 
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Hence,  the  first  roof  line  integral  can  be  evaluated  as 


K1  =  y,1(x',y')\H\Mx',y')) 

=  h?  —  [  [  e~^  -^5(y') -2  sm9^(9(y')  -  9(-y'))  sec Odx'dy' 

2m  Jx,=b  Jy,=0_e  V  b  )  [  b  b  J 

—h2  A2  —2a  J  .  f+e  (  |y'|\  T-2  „  2 a  sin 9  ,s  ,  , 

/  e  b  dx  / 


1-M) 


-h2H2  r  & 


2m  2  a 


C,  (-¥)  >»  -  “  r  (>  -  ¥)  <«  -  H  " 


-h2  A2  sec  9  f  b 


2m  [2  a 
K2  A 2  e~2a  sec  9 


r2ai  I--  -  —  (0) 

11  b  V  y  ’ 


(178) 


For  the  left  arm,  x'  =  —  oo  to  x'  =  —  b,  y'  —  0  and  the  trial  function  near  its  proximity  is 
given  by 

^2  =  a(  (179) 


'h  =  A\1~b) 

Using,  \y'\=y'  [%')  -  9{-y%  and  %')  =  S(-y'),  we  get 

=  TW) 

A  b  _  At  e¥  =  “  (i  _  At  e¥ 

dx'  \  b  )  b  I  b  J 


dx'dy ’ 


fT  ,  A  -b2  r-2  C •/  A  2a  SUld 

h^2  =  a^  -ys(y)  +  — p — (% ) - 0(-^/ ))  e  6 


Hence,  the  second  roof  line  integral  can  be  evaluated  as 


I<2  =  (4’2(x',y')\H\Mx',y')) 

=  zlid!  r  r  (i  _  Mi 

2m  J x/=_00  Jyi—Q_e  y  b  J 
X  zl5(y')  +  ^^  (e(y')-e(-y'))  sec  9  dx’ dy’ 
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- h2A 2  rb 

2  m  .L,-_ 


2a  x 1  . 

e  b  dx 


(lO-x 


x  -2  6( y')  +  (»(</)  -  #(-»'))]  dy' 


(180) 


Using  transformation  u  =  —x',  so  that  du  =  —dx' ,  therefore,  when  x'  =  — b ,  u  =  b  and  when 
a;'  =  —  oo,  u  =  oo. 

—ft2 A2  sec  6  fb  -2cu  f+e  f  It/'I  \ 

Aa= — ^Le^duL  A1 

-2..  .  2asind  ..  /1/  ,  . 

x  ~y%)  + — p — (%)  -^(-1/))  dV 


=  -^sec9  re^dum  0_M 
2m  .yu=6  A/=-e  V  o 

X  ±W)  +  (»(</)  -  »(-!/))  A/ 

—h2  A2  sec  6  T  6  oJ 


r»oo  /*+€ 

—  2au  I 

ebdu 


2m  2  a 


It/'I  \  —2  ,  ,  2a  sin  6  f+e 

^  T W-V  +  -P-L 


-h2  A2  sec  6  b  _2a  2  2ck 

2m  2a6  b  62 


1-y  (%')-W) 


h2  A2 e  2a  sec  0 


= 


(181) 


For  the  top  arm,  y'  —  b  to  y'  —  oo,  x'  —  0  and  the  trial  function  is  given  by 


(182) 


As  there  is  no  discontinuity  in  the  derivative  of  ^3  with  respect  to  y' ,  contribution  is 
not  taken  in  the  calculation  of  matrix  element  of  H  in  this  roof  line  as  it  is  already  taken 
into  account  in  the  bulk  integration.  Hence,  we  have  to  consider  H  as 


_  -h2  (  d2  ^  .  n  d2  \ 

2m  \dx'2  dx'dy' ) 

For  the  bottom  arm,  1/  =  —00  to  y'  =  —b,  x'  —  0  and  the  trial  function  is  given  by 

V>3  =  A  I  1  ~  —  J  e  * 


(183) 


(184) 
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With  simple  analogy  with  K\  and  K2  integrals,  it  is  clear  that  K:i 
Hence, 

4 

kj  =  y,k* 

i=l 

=  4  K  i 

_  2  h2  A2  e~2a  sec  9 
m  a 


I\\  and  K4  =  K\ . 


(185) 


B.  Calculation  of  matrix  elements  for  inner  four  arms 

The  inner  four  arms  are  defined  by  the  following  equations: 

x'  =  0  to  x'  —  b,  y'  —  0 

x'  =  —b  to  x'  —  0,  y'  —  0 

y'  =  0  to  y'  =  b,x'  =  0 

y  =  —b  to  y'  —  Qjx'  —  Q 

For  the  inner  right  arm,  x'  =  0  to  x'  =  b,  y'  =  0  and  the  trial  function  is  given  by 

-  A  ^1  -  e~a  (186) 

As  there  is  no  discontinuity  in  the  derivative  of  Ti  with  respect  to  x',  contribution  is 
not  taken  in  the  calculation  of  matrix  element  of  H  in  this  roof  line  as  it  is  already  taken 
into  account  in  the  bulk  integration.  Hence,  we  have  to  consider  H  as 

_^2  /  Q2  Q2  \ 

H  =  2^T  \W2  ~  2  Sm  6  dx’dy’ )  ^187') 

Note :  Theta  functions  are  the  special  class  of  functions  which  is  defined  by 

9(x  —  x')  —  1  when  {x  —  x')  >  0 
=  0  when  {x  —  x')  <  0 

Also,  another  important  property  of  theta  function  is  it’s  derivative  is  Dirac-delta  function, 
that  is 


—  (9(x  —  x'))  =  5(x  —  x') 
ax 
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and  Dirac-delta  function  is  defined  by 

8(x  —  x')  —  1  when  x  —  x'  =  0 
=  0  when  x  —  x'  ^  0 

In  the  following  calculations,  we  make  use  of  these  properties  of  both  functions  to  arrive  at 
the  result. 


Using,  | y'\=y'  [6{y')  -  0(-y%  and  8{y')  =  8(-y'),  we  get 


<92^i 


<92^i  2  Ax'. 

^T  =  -— %)e 

(9^1  Ae~a  , 


=  A- 


Ae  a 

“d  ^ =  — 

U  feV)  +  ^  W)  -  «(-</'))]  e- 

\m  b  b- 


Hence,  for  this  case  the  roof  line  integral  can  be  evaluated  as 


L\  =  (^>i(x'  ,y')\H\^i(x'  ,y')) 

-k2A2e~2a  fb  f0+e 


Jx'= 0  Jy'= 0-e  V  u  / 

— 2  sin  $ 

■io-5(l/)  +  (%')  -  «(-!/))  sec^W 

£r  o-  ' - „ - / 


-h2A2e  2"sec 9  [b  \—2x'  2  sin# 


—  h2A2e~2a  sec  6 

2m 

—  h2A2e~2a  sec  6 

2m 

h2A2e~2a  sec  6 


(0)  dx 


—  /  x'dx'  +  0 

b  )  Jrrl—Ci 


seed  r /  — 2 \  b 2 
1 1?)  ~2 


(188) 


For  the  inner  left  arm,  x'  =  —b  to  x'  =  0,  y'  =  0  and  the  trial  function  is  given  by 


^2  —  A  \  1  + 


(189) 


Using,  | y'\=y'  \6{y')  -  0(-y%  and  8{y')  =  8(-y'),  we  get 
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<92T2  2  Ax' 


5{y')e 


d^2  _  Ae~a 
3x'  b 2  ^ 


<92T2  _  Ae~a 
dx'dy'  b2 


W)  -  «(-</)] 


-h2  \2x' 


H^  =  A^[irs^ 

The  roof  line  integral  can  be  obtained  as 


2  sin  6 


0{-y'))  e 


L2  =  ^2(x',y')\H\^2(x/,y')) 

-h2A2e~2a  r°  r0+e 


Jx'=-b  Jy'= 0-e 
2x'  f/.A  2sil10  //./.A 


x  P%')  -  PP  («(!/)  -  PP  sec  6 dx'dy1 
tr  bz  <« - v - ' 

— h2A2e~2a  sec  0  /‘°  2x'  2  sin  0 

2m  .L'-h  b2  b2 


(0)  da;' 


Using  substitution  u  =  —x',  so  that  du  =  —da;',  therefore,  when  x'  =  — 6,  w 
a;'  =  0,  u  =  0. 


=  b  and  when 


h2  A2  e  2a  sec  6 


h2A2e~2a  sec  0  /  2  \  fb 


u  du  +  0 


2m  [\b2J  Ju=0 

h2A2e~2a  sec  6  f/2  \  b2' 

2  m  [\j?)  ~2  _ 

h2A2e~2a  sec  9 


(190) 


For  the  inner  top  arm,  y'  —  0  to  y'  —  b,  x'  —  0  and  the  trial  function  is  given  by 


#3  =  A  1 


PA  -a 

b2  ) 


(191) 


In  the  same  analogy  as  integral  L i,  it  can  be  shown  that  integral  for  the  inner  top  arm 
L3  =  L\.  For  the  inner  bottom  arm,  1/  =  —6  to  y'  =  0, a;7  =  0  and  the  trial  function  is  given 
by 


+  e-« 


(192) 
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With  the  same  calculations  as  integral  L2,  it  can  be  shown  that  integral  for  the  inner  bottom 
arm  L4  =  L2.  Hence, 

4 

Kn  =  ^  Li 

i=  1 

=  4  Lx 

_  2  h2  A2  e~2a  sec  6 
m 


C.  Calculation  of  matrix  elements  for  outer  eight  sides 

The  outer  eight  arms  are  mathematically  defined  by  the  following  sets  of  equations: 

y'  =  0  to  y'  =  b,  at  x'  =  b 

y'  =  —b  to  y'  =  0,  at  x'  =  b 

y  =  0  to  y  =  b,  at  x'  =  —b 

y'  =  —b  to  y'  =  0,  at  x'  =  —b 

x1  =  0  to  x'  =  b,  at  y  =  b 

x1  =  —b  to  x'  =  0,  at  y  =  b 

x'  =  0  to  x'  =  b,  at  y'  =  —b 

x1  =  —b  to  x'  =  0,  at  y  =  —b 

For  the  side,  y'  —  0  to  y'  —  b  at  x'  =  b,  the  trial  function  is  given  by 

(pa  =  A^l-  e~a  for  x'  <  b 

=  A  ^1  —  for  x'  >  b 

The  Hamiltonian  operator  is  given  by 

—  fjf  /  52  Q2  Q2  \ 

H=^  [d^  +  d^~2  sm9dx'dy' )  ^194'1 

As  there  is  no  discontinuity  in  the  derivative  of  ipa  with  respect  to  y' ,  contribution  from 
is  not  taken  in  the  calculation  of  matrix  element  of  H  for  this  side  as  it  is  already  taken  into 
account  in  the  bulk  integration.  If  this  term  in  Hamiltonian  operator  is  not  dropped  then 
we  get  the  extra  contribution  which  is  incorrect.  Hence,  we  have  to  consider  H  without  the 
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second  term  as 


H=  —  (--2  ■  &2 

2m  ^cU'2  Sm  dx'dy' 


Using  ^-function,  the  trial  function  can  be  written  as 


<Pa  =  A  (l  —  J  9(x'  -b)  +  j6(b-  X ')  j  e 


-x')+f  e(x’-b) 


(195) 


(196) 


^  =  A\  5(x'  -b)-  ^ 5(6  -  x')  +\e(b  -  x') 

r)  yU  h  A  A 


+  a(i-^  9(x'  -b)  +  jd(b  -  x')  )  e-«[«(6-*')+^®(*'-6)] 

CV  CV 

x  aS(b-x')  -  —  5(x' -b)--9(x' -b) 
b  b 

" - V - ' 

L  =0  J 

=  ~  \yT9(b  -  x')  +  aQ{x'  -  6)  (l  ~  xll  e~a\.0ib~x,)+^0{x'~b)'\  (197) 

b  l b  \bj 

g(x',  y')  =  j9{b-  x')  +  a  9(x'  -b)  ^1  - 
and  f(x')  =  —a  9(b  —  x')  +  y  9{x'  —  b) 

K  =  _i>S(b-,)+aS(x._i)U_l\ 


f)  f  Tf  I 

^-  =  -a  s(b  -  X*)  +  Six'  -b)+y9(x'  -  b) 
ox  .  0  .0 


Aef{x,)  S(b  —  x') 
0 


—  CM 

=  —  9(x'  -  b ) 
b 

Therefore  9f'!  =  —  ^ g(x' ,  y')  e^x'^ 
dx'  b  ,yy 

d^  =  _A\dg_  df\  w 
dx'2  b  dx'  y  dx' 


y'  ,  /-■  y\]  «2  n  y\Q{  >  U\ 

~  +  a(i~j)\  -6) 

lntegral=0  / 

=  --e~a)5(x' -b)  a  —  (a  +  1)  — 
b  b 


(198) 


Also, 


=  -±JL  \q(x'  „')e/<«') 

dx'dy'  b  dy'  l 


=  _Aef{x,)d1_ 


=  -^e/(a/)  [^0(6  -  x')  -  X>  0(x'  -  b ) 
b  [b  b 


=  -^e/(a;,)  9(b-x')-a9(x' -b) 
bz  v - v ' 


lntegral=0 


Therefore 


—li2  [A  /  7/ 

—  — -e-^fx' -  6)  (a- (a +  1)4- 

2  m  6  v  ’  V  v  ;  6 


(199) 


Hence,  for  the  side  at  x'  =  b  extending  from  1/  —  0  to  1/  —  b,  the  integral  is  given  by 


Si  =  {(Pa(x,,y,)\H\(pa(x,,y1)) 

fc2  fb  rb+e  f  A  / 

= -—  /  /  <Pa(x',y')  —r  e~a5{x'  -  b)  (a  -  (a  +  1)  y- 

2  m  .7)/=0  .7x/=6-e  L  6  V  b 


h2  sec  9  (  A ' 


1  r 


h2  sec  9  (  A ' 


b  r 


<Pa(b,y')  a  —  (a  +  1)  —  dy' 

a(  1  —  e~a  a  —  (a  +  1)  ^  dy' 


h2A2e  2a  sec  9  fb  (  y'\  \  y'~\  , 

2 mb  L  V-j)  r-(a+i)ij dy 

2^  [a  -  (2a  +  i)  -  +  (a  +  i)  — J  dy 

h2  A2  e~2a  sec  9  \  b2  b3 

—  |a6-(2a:+l)-  +  (a+l)  — 

h2A2e~2asec9  .  .b  .  .6] 

—  ^6-(2a+l)-  +  (a  +  l)-^ 

/i2  A2  e_2Q  sec  9  \  2a  +  1  a  +  1 

-  OL  —  -  - 

2m  [  2  3  J 

h2  A2  e~2a  sec  9  fa  1\ 


2  m  \ 3  6  j 

h2  A2  e~2a  sec  9  (2  a  —  1 


sec  9  dx'  dy' 


(200) 


Use  of  appropriate  trial  functions  for  the  other  sides,  it  can  be  shown  that  the  integrals  Sn, 
where  n  =  (1,  8^  are  equal.  Hence,  the  total  contribution  to  the  integrals  from  the  eight 
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sides  is  obtained  as 


Km  =  8  Si 

h 2  A2  e~2a  seed  / 2  a  —  1 

2  m  \  6 

2h2  A2e~2asec9 
=  - o -  (2  a  -1) 


(201) 


The  matrix  element  of  H  for  all  the  roof  lines  is  given  by 


(H{a))Rooflinea  =  KI  +  KII  +  KIII 

2h2  A2e~2asec9  2  h2  A2  e~2a  sec  9  2  h2  A2  e~2a  sec  9  , 

=  - + - + - - -  2a 

ma  m  ■>  m 

_  2  h2  A2  e~2a  sec  9  /12a-l\ 
m  \  a  3  / 

2  h2  A2  e~2a  sec  9  f  2  a2  +  2  a  +  3\ 
m  \  3  a  J 


1) 

(202) 


